AD-A055  857  FLORIDA  STATE  UNIV  TALLAHASSEE  GEOPHYSICAL  FLUID  DYN~ETC  F/G  4/2 
CLIMATE  SIMULATION. (U) 

MAR  78  R L PFEFFER  N00014-76-C-0844 

UNCLASSIFIED  NL 


ADA055857 


Principal  Investigator 
'^Geophysical  Fluid  Dynamics  Institute 
Florida  State  University 


18  06  IS  100 


‘-/t! 


DISTRIBUTION  STATEMENT  K 

Approved  for  public  teleosel 
Distribution  Unlimited 


c. 


i 


f 

i 


PREFACE 

Under  contract  N00014-76-C-0844  (NR  062-546)  field  distributions  of 
fluid  temperature  and  speed  were  measured  in  thermally-driven  rotating 
fluids  which  simulate  important  aspects  of  meteorological  and  climatic 
atmospheric  behavior.  Climatic  statistics  were  generated  which  have 
relevance  to  problem  of  extended  range  atmospheric  prediction.  The 
results  of  our  experimental  work  and  their  interpretation  are  contained 
in  two  published  papers  which  comprise  this  Final  Report. 


A.  Cyclic  Variations  . of  the  Imposed  Temperature  Contrast  in  a 
Thermally  Driven  Rotating  Annulus  of  Fluid,  J.  Atmos.  Sci.,  35,  1978  by 
Buzyna,  Pfeffer  and  Kung;  This  paper  presents  the  results  of  an  extensive 
series  of  experiments  with  periodic  "seasonal"  temperature  changes  super- 
imposed upon  otherwise  known  flows  in  the  range  of  parameters  pertinent 

i 

to  atmospheric  weather  and  climatic  variations.  In  these  experiments  the 
rotation  rate  (R)  was  held  constant  and  the  imposed  temperature  difference 
across  the  fluid  (AT)  was  varied  in  a continuous,  periodic  fashion  for  up 
to  10  "years"  after  an  initially  stable  climate  was  established.  Nhereas 
the  flow  in  conventional  annulus  experiments  was  found  to  be  deteznCna.t» 
in  the  sense  that  it  depended  only  upon  the  point  in  dimensionless>para- 
meter  space  at  which  the  experiment  was  conducted  and  the  direction  frost 
which  it  was  approached,  the  flow  configuration  in  the  experiments  in  which. 
AT  was  varied  cyclically  was  found  to  be  indeterminate  in  the  sense  that  it 
changed  from  one  state  to  another  sporadically  and  could  not  have  been  pre- 
dicted solely  as  a function  of  "season".  During  different  annual  cycles 
the  general  circulation  of  the  fluid  in  a particular  season  took  on  dif- 
ferent, behavioral  characteristics.  For  ••xample,  the  flow  pattern  which 
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developed  in  one  summer  season  sometimes  established  itself  and  shaped  events 

in  the  following  winter  season.  Sometimes  the  patterns  persisted  into  a 

third  year  before  breaking  into  a new  regime  (e.g.  wave  number  change,  or 

change  from  amplitude  vacillation  to  wave  dispersion).  This  is  similar  to 

atmospheric  events  in  which  synoptic  waves  become  established  in  different 

phases  from  normal,  or  in  which  the  jet  stream  becomes  established  north  - 

% 

or  south  of  its  normal  position  for  considerable  periods  of  time.  Such 

/ 

abnormalities  in  the  atmosphere  are  often  accompanied  by  extended  periods 
of  drought,  altered  hurricane  paths  and/or  early  onset  of  trinter  conditions 
such  as  were  experienced  in  recent  winters.  There  are,  of  course,  a number 
of  feedback  mechanisms  at  work  in  the  atmosphere  (including  cloud  and  ice 
reflectivity,  atmosphere>ocean  interactions,  mountain  influences,  etc). 

The  fact  that  important  similarities  exist  between  atmospheric  behavior  and 

i 

the  behavior  of  fluids  in  rotating  annulus  experiments  in  which  ocean,  land 
and  cloud  irregularities  are  not  present,  however,  emphasizes  the  fundamental 
role  of  the  annulus  model  in  studies  of  non-linear  processes  which  shape 
climate  and  climatic  changes. 

B.  Catastrophic  Changes  in  Circulation  Flow  Patterns,  American  Scientist, 
65  (5),  614-621,  1977  by  Lacher,  McArthur  and  Buzyna:  One  by-product  of  our 
"seasonal"  experiments  was  the  observation  that  strong  hysteresis  effects 
are  present  when  the  imposed  AT  is  varied  very  slowly.  In  particular, 
certain  wave  number  changes  take  place  at  different  points  in  dimensionless- 
parameter  space  depending  upon  whether  AT  is  slowly  increased  or  decreased. 
This  phenomenon  is  one  which  was  reported  on  earlier  by  Fultz  and  others. 
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In  our  laboratory  we  have  conducted  a series  of  such  experinents  with  multi- 
probe temperature  measurements  and  have  determined  the  available  potential 
energies  associated  with  the  two  states  possible  at  each  of  several  points 
in  dimensionless-parameter  space.  The  important  dimensionless  parameters 
for  these,  flows  are  the  thermal  Rossby  number  (ratio  of  inertial  to  Coriolis 
acceleration),  the  Taylor  number  (square  of  the  ratio  of  the. Coriolis  accele- 
ration to  the  visco\is  force)  and  the  Prandtl  number  (ratio  of  the  viscous 
to  the  thermal  dissipation  coefficient) . We  found  that  dependi^  tqion  the 
direction  from  .which  a point  in  dimensionless-parameter  space  was  approadxe<» 
different  energy  levels  exist.  This  led  to  the  paper  by  Lacher,  McArdtur 
and  Buzyna  in  which  they  presented  an  alternative  method  of  representi^ 
bimodal  hysteresis  pheonomena  using  Thom's  catastrophe  theory  (see  his  book 
entitled  "Structural  Stability  and  Morphogenesis"  Reading  MA:  W.  A.  Benjamis 
Inc.  1975).  The  model  proposed  is  the  simplest  possible  TOdel  for  describ- 
ing hysteresis  and  intransitivity  phenomena  observed  in  laboratory  experiments 
displaying  abrupt  climatic  changes.  The  behavior  of  the  fluid  in  such  ex- 
periments is  foimd  to  be  closely  approximated  by  a cusp  model  of  the  form 

t = f(fl,  AT,  E) 

where  E is  a measure  of  the  available  energy  in  the  system,  £ is  the  time 
rate  of  change  of  E,  n is  the  rate  of  rotation  of  the  system  and  AT  is  the 
imposed  "pole  to  equator"  temperature  difference.  The  function  f is  the 
)(vcdlent  of  a potential  function  whose  qualitative  properties  are  not 
changed  by  small  perturbations.  A reprint  of  this  paper  with  further 
details  is  contained  at  the  end  of  this  Final  Report. 
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1.  Introduction 

The  earth's  atmosphere  is  subjected  to  periodic  thermal 
forcing  with  a cycle  period  of  365  1/4  rotations.  While  varia- 
tions of  the  zonal ly- averaged  temperatures  in  the  tropics  are 
small,  arctic  and  mid-latitude  temperatures  undergo  marked 
seasonal  oscillations.  In  winter,  when  arctic  temperatures 
are  lowest,  intense  baroclinic  disturbances  are  in  evidence 
over  the  hemisphere  and  the  hemispheric  average  eddy  available 
potential  and  kinetic  energies  reach  their  peaks.  In  summer, 
the  eddy  energies  reach  their  lowest  magnitudes.  It  is  not 
obvious,  a priori,  whether  more  intense  synoptic  wave  activity 
is  universally  associated  with  greater  thermal  forcing,  or 
whether  this  result  is  a peculiarity  of  the  earth's  atmospheric 
circulation  due  to  the  particular  rate  of  rotation  of  this 
planet  and  the  range  of  intensities  of  the  thermal  forcing  to 
which  it  is  subjected.  Indeed,  laboratory  experiments  and 
baroclinic  theory  reveal  that  synoptic  waves  are  replaced  by 
symmetric  flow  when  the  imposed  thermal  forcing  exceeds  a 
certain  critical  magnitude. 

In  spite  of  the  periodic  nature  of  the  thermal  forcing 
function,  the  scale,  phase  and  propagation  speed  of  the  syn- 
optic disturbances  display  intra-annual  differences  which,  in 
turn,  are  accompanied  by  rather  signific2mt  differences  in 
weather  phenomena  in  the  same  season  of  different  years.  It 
is  recognized,  in  this  connection,  that  the  earth's  atmosphere 
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is  also  subjected  to  geographically  irregular  thermal  and  oro- 
graphic influences  and  that  it  participates  in  a complex  set 
of  feedback  mechanisms  involving  releases  of  latent  heat  of 
condensation  and  sublimation,  exchanges  of  heat  with  the  oceans, 
cloud  and  ice  reflectivity,  etc.  At  present,  it  is  not  yet 
known  to  what  extent  intra-annual  differences  in  weather  patterns 
are  attributable  to  the  complexities  of  the  forcing  functions 
and  feedback  mechanisms,  and  to  what  extent  they  are  inherent 
in  the  system  due  to  the  non-linear  property  of  hysteresis 
("memory" ) . 

It  is  the  purpose  of  the  present  article  to  explore  answers 
to  such  questions  using  a laboratory  experimental  model  of  a 
rotating  fluid  subjected  to  periodic  thermal  forcing.  The  di- 
mensionless parameters  of  the  model  are  controlled  such  as 
to  produce  quasi-geostrophic,  quasi-non-divergent , strongly 
stratified  flow  similar  to  that  in  the  earth's  atmosphere.  The 
model  is,  however,  by  design, devoid  of  many  of  the  complexities 
of  the  atmospheric  forcing  function  and  of  the  feedback  mechan- 
isms to  which  the  atmosphere  is  subjected. 

Snyder  and  Youtz  (1969)  have  made  a limited  study  of  the 
behavior  of  such  a fluid  subjected  to  a sinusoidally  varying 
Imposed  temperature  contrast  in  a region  of  dimensionless- 
parameter  space  centered  on  the  static  transition  between  wave 
numbers  2 and  3.  The  region  in  which  their  experiments  were 
conducted  displayed  weak  hysteresis  properties.  With  roughly 
a 30%  variation  around  a nominally  10®C  imposed  temperature 
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contrast,  they  found  that,  in  the  limit  of  slow  cycling,  transi- 
tions between  wave  numbers  2 and  3 tended  to  take  place  near  the 
mean  value  of  the  imposed  temperature  contrast  for  such  transi- 
tions under  static  conditions.  With  decreasing  cycle  period 
the  transition  points  deviated  increasingly  from  the  static 
value  until,  at  sufficiently  small  cycle  periods,  no  transition 
took  place. 

In  the  present  study  experiments  were  conducted  at  three 
different  rotation  rates,  with  four  different  ranges  of  varia- 
tion of  the  imposed  temperature  contrast,  and  with  four  different 
cycle  periods.  Temperature  spectra  and  eddy  and  zonal  tempera- 
ture variances  were  calculated  and  are  presented  in  subsequent 
sections  as  a function  of  time.  The  region  of  dimensionless- 
parameter  space  in  which  these  experiments  were  conducted 
displays  strong  hysteresis  properties. 

In  order  to  provide  a proper  framework  for  interpretation 
of  the  results  of  experiments  with  periodically  imposed  thermal 
forcing,  we  shall  review  first  some  of  the  properties  of  thermal- 
ly driven  rotating  flows  which  develop  when  the  thermal  driving 
force  is  not  varied  with  time.  We  shall  also  examine  the 
behavior  of  fluids  subjected  to  small  step  changes  of  the  thermal 
forcing  function.  Attention  will  be  focused,  in  particular,  on 
(1)  the  levels  of  thermal  variance  attained  by  different  synoptic 
wave  components  at  relevant  points  in  dimensionless-parcuneter 
space,  (2)  the  response  time  of  the  fluid  to  step  changes  in  the 
thermal  forcing  function,  and  (3)  the  hysteresis  characteristics 
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displayed  when  a point  in  dimensionless-parauneter  space  is 
approached  from  different  directions  (i.e.,  by  lowering  or 
raising  the  thermal  forcing  function) . 


2.  List  of  symbols  eind  dimensionless  parameters 


R radial  distance  from  the  cucis  of  rotation  to 
the  annulus 

z height  above  the  base  of  the  fluid 
t time 

a radius  of  inner  cylinder  (7.303  cm) 

b radius  of  outer  cylinder  (14.923  cm) 

depth  of  working  fluid  (15.0  cm) 

T azimuthal  ("zonal")  average  of  the  fluid  temp- 
erature at  a specified  radius  ("latitude") 
inner  bath  temperat\ire 

AT  imposed  temperature  contrast  between  the  inner 
and  outer  cylinders  containing  the  working 
fluid 

rate  of  rotation  of  annulus 
g acceleration  of  gravity 

a coefficient  of  volume  expansion  of  the  fluid 
p fluid  density 

V coefficient  of  kinematic  viscosity  of  the  fluid 
k coefficient  of  thermometric  diffusivity  of  the 
fluid 

Pr=y/k  Prandtl  number  (57)  2 2 

Rom  imposed  thermal  Rossby  number  {=gaHLT/lQ  (b-a)  ]} 

2 i»  2 

Ta  Taylor  number  [=4(2  (b-a)  /v  ] 

VgfVg  zonal  and  eddy  temperature  variance,  respectively, 
at  mid-depth  in  the  fluid 
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3.  Apparatus  and  procedures  ; 

I 

The  working  fluid  used  in  this  experimental  series,  5 
centistoke  silicone  fluid,  depth  15  + 0.05  cm,  was  contained 
in  an  annular  region,  the  inner  and  outer  radii  of  which  were 
7.303  and  14.923  ± 0.005  cm,  respectively.  The  annular  region 

J 

was  formed  by  two  concentric  brass  cylinders  of  0.43  cm  wall  j 

thickness  and  20.3  cm  length  imbedded  in,  and  oriented  normal 
to,  a flat  Plexiglas  base  2.4  cm  thick.  The  cylindrical  walls 
were  maintained  at  prescribed  temperatures  by  means  of  circu- 
lating water  baths  in  contact  with  the  outside  of  the  outer 
wall  and  the  inside  of  the  inner  wall.  The  water  in  each  bath 
was  circulated  at  a rate  of  2 gal.  min~^  and  its  temperature 
was  controlled  to  within  ± 0.05®C. 

The  entire  annulus  assembly  was  mounted  on  a turntable  with 
a hydrostatic  oil  bearing  with  its  axis  of  symmetry  concentric 
with  the  vertically-oriented  rotational  axis  of  the  turntable. 

The  turntable  was  driven  by  a direct-drive  servo-controlled  dc 
torque  motor  which  provided  long-  and  short-term  rotational 
stability  of  the  order  of  ± 0.01%.  The  axis  of  rotation  was 
aligned  within  5 sec  of  arc  with  the  vertical  (gravitational 
field) . 

The  internal  distribution  of  fluid  temperature  was  measured 
with  an  array  of  bead  thermistors  positioned  within  the  fluid  by 
suspending  them  by  their  long  leads  (diam:  0.005  cm).  The 
entire  observational  network  embedded  in  the  fluid  consisted  of 
a distribution  of  14  single  thermistors  for  temperature  measure- 





merits  and  21  "quadruple"  velocity-temperature  transducers  located 
at  different  depths  and  radii  in  the  fluid  as  shovm  in  Fig.  1. 
Each  quadruple  consisted  of  three  bead  thermistors  forming  an 
equilateral  triangle  in  the  horizontal  plane  which  sensed  the 
temperature  field  produced  by  a fourth  heater  bead  located  at 
the  center  of  the  triangle  (as  described  by  Fowlis,  Buckley  and 
Ruppert,  1972) . These  were  present  in  the  fluid  for  the  purpose 
of  calibration  and  testing.  Only  temperature  measurements  were 
used  in  the  present  study.  The  total  number  of  wires  passing 
vertically  in  the  fluid  (including  return  leads)  was  154.  The 
details  of  the  thermistor  design  and  circuity,  the  digital  data 
acquisition  system  and  the  procedures  used  to  extract  and  process 
the  data  are  described  by  Pfeffer,  Buzyna  and  Fowlis  (1974). 

The  cyclic  variation  of  the  temperature  contrast  (AT)  across 
the  annulus  walls  was  accomplished  by  cycling  the  cold  inner-bath 
temperature  (T^ ) while  the  warm  outer  bath  temperature  was  main- 
tained constant.  A modified  slider  crank  mechanism  was  used  to 
produce  a cyclic  variation  of  the  thermoregulator  which  control- 
led the  cold  bath  temperature.  In  this  version 

CT  the  connecting  link,  of  fixed  length,  was  a string 
which  wound  around  a shaft  connected  directly  to  the  thermo- 
regulator (mecury  thermometer  type  with  a rotating  head) . The 
turning  of  the  thermoregulator  varied  the  cold  reservoir  bath 
temperature  in  a cyclic  fashion.  ' The  ratio  of  the  distance 
between  the  thermoregulator  shaft  and  the  crank  pivot  to  the 
length  of  the  crank  was  5.5.  This  combination  gives  a maximum 


Fig.  1: 


Sketch  showing  horizontal  distributions  of  thermistor  probes  at 
two  levels  in  the  annulus.  The  concentric  circles  represent  radii 
R = 8.57,  11.11  and  13.65  cm  (1/6,  1/2  and  5/6  gap-width, 
respectively).  The  locations  of  single  thermistors  at  mid-depth 
(7.50  cm  from  the  base  of  the  annulus)  are  designated  by  solid 
circles  (•) . The  locations  of  "quadruples"  (transducers  consisting 
of  four  thermistors  used  for  measuring  fluid  temperature  and  velocity) 
at  7.50  and  11.25  cm  above  •’he  base  of  the  annulus  are  designated 
by  solid  (a)  and  open  (a)  triangles,  respectively. 
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deviation  from  sinusoidal  variation  of  approximately  4.5%  at 
points  mid-way  between  the  extreme  temperatures.  Furthermore, 
the  shape  of  the  cyclic  curve  is  broader  than  sinusoidal  in  the 
region  of  the  warm  end  of  the  bath  temperature  and  narrower  at 
the  cold  end.  Hence  the  duration  of  the  minimum  temperature 
contrast  is  slightly  longer  than  the  maximum  temperature  contrast. 

4 . Experiments  conducted 

The  region  in  dimensionless-parameter  space  within  which 

the  present  experiments  were  conducted  is  shown  in  Fig.  2.  The 

dots  and  the  transition  curve  separating  the  symmetric  regime 

from  the  wave  regime  correspond  to  a previously  unpublished 
2 

series  of  null  experiments  using  a two-level,  100-probe  network 
of  temperature  measurements.  Details  of  the  network  and  exper- 
imental procedures  are  described  by  Buzyna,  Kung,  Fowlis  and 
Pfeffer  (1973) . The  numbers  next  to  each  dot  in  the  figure 
designate  the  dominant  wave  forms  (in  order  of  descending 
amplitude)  observed  in  each  experiment.  Nearly  all  of  the  wave 
experiments  exhibited  some  degree  of  amplitude  vacillation  of  a 
single  dominant  mode,  or  of  three  adjacent  dominant  modes  with 
different  mean  eimplitudes  and  different  magnitudes  of  vacillation. 
The  single  exception  was  an  experiment  designated  in  Fig.  2 by 
4,  5,  6D  which  exhibited  wave  dispersion  without  significant 
variations  in  the  amplitudes  of  the  dominant  waves.  Three  addi- 
tional null  experiments,  indicated  by  crosses  in  Fig.  2,  were 


2 

Null  experiments  are  defined  here  as  experiments  with  the 
externally  controlled  parameters  (rotation  rate  and  imposed  tem- 
perature contrast)  held  fixed. 
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Fig,  2:  Regime  diagram  showing  the  dominant  wave  numbers  (in  order  of 
descending  amplitude)  which  form  at  different  values  of  the 
thermal  Rossby  number  (Roj)  and  Taylor  number  (Ta)  in  an  annulus 
of  silicone  fluid  (Prandtl  number  57)  with  depth  H = 15.0  cm  and 
radial  dimensions  7.30  ^ R 14.92  cm  when  the  imposed  temperature 
contrast  (AT)  and  the  rotation  rate  (fi)  are  held  constant.  The 
data  next  to  the  dots  are  previously  unpublished  results  obtained 
with  the  same  fluid  in  the  same  annulus  as  in  the  present 
experimental  series,  with  a different  distribution  of  100 
thermistors  at  two  depths  in  the  fluid.  The  two  entries  at 
RO.J.  = 0.50,  Ta  = 1.43  x 10^  indicate  the  two  different  wave 
numbers  which  formed  when  the  same  experiment  was  run  on  different 
days.  The  data  next  to  the  crosses  are  the  results  of  additional 
null  experiments  conducted  as  part  of  the  present  investigation 
with  the  array  of  probes  shown  in  Fig.  1.  The  vertical  line 
segments  labeled  A,  B,  C,  D and  E denote  the  ■ ranges  over 

which  the  imposed  temperature  contrast  was  cycled  in  the  present 
experimental  investigation  (see  Table  1) . 
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conducted  with  the  present  observational  network  to  help  clarify 
the  nature  of  the  different  wave  states  encountered  in  experi- 
ments with  cyclic  variations  of  the  imposed  temperature  contrast. 
The  control  parameters  and  the  associated  dimensionless-parameters 
corresponding  to  these  experiments  are  listed  as  experimental 
series  G in  Table  1. 

Fig.  3 shows  the  magnitude  of  the  time-mean  eddy  thermal 
variance  of  the  fluid  at  the  different  points  in  dimensionless- 
parameter  space  at  which  null  experiments  with  the  earlier  100- 
probe  network  were  conducted.  Since  imposed  periodic  fluctua- 
tions of  the  temperature  contrast  (AT)  represent  vertical  paths 
in  the  Ro^  vs.  Ta  plane,  it  is  clear  by  inspection  of  this 
diagram  that  whether  or  not  the  wave  amplitudes  will  be  larger 
in  "winter"  (high  imposed  temperature  contrast)  than  in  "summer" 
(low  imposed  temperature  contrast)  will  depend  upon  where  in 
dimensionless-parameter  space  the  experiment  is  conducted. 

Judging  by  the  observed  behavior  of  the  earth's  atmosphere,  it 
would  appear  that  experiments  C and  D (Fig.  2 and  Table  1)  should 
be  the  most  relevant  of  the  present  series,  in  this  respect,  to 
atmospheric  seasonal  changes. 

In  addition  to  the  null  experiments,  experiment  SE-32, 
series  F (Table  1)  was  conducted  to  determine  the  response  char- 
acteristics of  the  fluid  tc  a near-step  change  in  the  cold-bath 
temperature  (T^)  of  4.6®C  with  n held  constant.  In  two  other 
experiments  in  series  F,  SE-34  and  SE-35,  T.  was  varied  in 
smaller  discrete  steps  over  16-  and  24-hour  periods,  respectively. 
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Time-mean  eddy  thermal  variances  ^.corresponding  to  the  experiments 
denoted  by  dots  in  fig.  2.  The  variances  here  have  been  converted 
to  energy  units  (lO"^  ergs  cm' 2)  by  multiplying  by  the  coefficient 
used  in  defining  the  available  potential  energy  of  a liquid. 
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with  sufficient  time  allowed  after  each  step-change  for  the  fluid 
flow  to  reach  a non-transient  state  with  respect  to  the  imposed 
external  conditions.  The  purpose  of  these  experiments  was  to 
elucidate  the  hysteresis  characteristics  of  the  system.  Hystere- 
sis was  also  studied  by  means  of  experiment  SE-36,  series  F,  in 
which  was  varied  continuously  over  a 24-hour  period.  This 
variation  was  sufficiently  slow  to  permit  the  flow  to  keep  up 
with  the  resulting  changes  in  AT  without  displaying  significant 
transient  effects. 

Experiments  designated  by  the  letters  A,  B,  C,  D 
and  E in  Table  1 and  Fig.  2 were  conducted  with  the  rotation 
rate  held  fixed  and  with  more  rapid  cyclic  variations  of  the 
imposed  temperature  difference  produced  by  continuously  varying 
the  temperature  of  the  cold  inner  bath.  In  these  experiments 
several  different  fixed  rotation  rates,  and  several  different 
frequencies  and  amplitudes  of  the  imposed  temperature  difference 
were  employed. 

5.  Characteristic  time  scales 

In  order  to  facilitate  the  interpretation  and  understanding 
of  the  experimental  results,  we  consider  here  several  different 
characteristic  time  scales.  These  are  (a)  the  amplitude  vacilla- 
tion periods,  (b)  .the  periods  over  which  the  imposed  temperature 
contrast  was  varied,  (c)  the  response  times  of  the  fluid  (i.e. , 


the  e-folding  times  required  for  the  fluid  to  respond  to  step 
changes  in  T.^)  , and  (d)  the  lag  times  required  for  the  fluid  to 
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respond  to  cyclic  variations  in  T.. 

(a)  The  amplitude  vacillation  periods;  The  characteristics 
of  amplitude  vacillation  were  described  by  Pfeffer,  Buzyna  and 
Fowlis  (1974).  The  vacillation  period  is  an  internal  scale 
resulting  from  the  general  dynamics  of  the  flow.  For  given 
externally- imposed  parameters  (either  fixed  or  transient)  the 
vacillation  periods  of  successively  higher  wave  numbers  in  these 
experiments  were  successively  shorter  (see  Table  2).  This  was 
the  case  when  different  wave  numbers  coexisted  in  the  same  exper- 
iment, and  when  the  dominant  wave  number  changed  during  an  exper- 
iment . /iT  was  also  the  case  when  different  wave  numbers  fomed 


as  a consequence  of  repeating  an  experiment  at  the  same  point  in 


. r ; ; 

dimensionless-parameter  space.  For  excunple,  in  experiment  SE-29 


(series  A)  the  dominant  modes  occurring  at  different  times  with 
different  relative  magnitudes  were  wave  numbers  3,  4,  and  5. 

Their  characteristic  vacillation  periods  were  approximately  130, 


60,  and  53  rot,  respectively. ^ The  vacillation  period  of  each 
mode  was  also  observed  to  decrease  with  successively  higher 
rotation  rates.  For  excunple,  at  the  higher  rotation  rate  of 
experiment  SE-30  (series  E) , the  dominant  modes  were  again  wave 
numbers  3,  4,  and  5.  The  vacillation  periods  of  these  modes  in 
this  experiment  were  approximately  120,  50,  and  44  rot,  respec- 
tively. The  greatest  uncertainty  in  these  determinations  is  in 
the  vacillation  period  of  wave  number  3 because  this  mode  occurs 
with  sufficiently  large  amplitude  only  over  one  or  two  major 
amplitude  vacillation  cycles.  Hence,  the  vacillation  periods 
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reported  for  this  wave  nximber  should  be  regarded  only  as  crude 
estimates. 

(b)  The  periods  of  the  imposed  temperature  contrast ; 

These  periods  were  deliberately  set  to  be  long  in  comparison 

with  the  thermal  lag  time  of  the  fluid  and  with  the  vacillation 

periods  observed  in  null  experiments.  The  imposed  periods 

employed  were  1/2,  1,  2,  16,  and  24  hours.  The  16-  and  24-hour 
for  the  purpose  of  studying 

periods  were  selectedAhysteresis.  The  1/2  to  2 hour  imposed 
cycles  represent  periods  ranging  from  492  to  2516  rot,  depending 
upon  the  rotation  rate  in  each  experiment  (see  Table  2) . In  most 
experiments  these  periods  were  more  than  ten  times  longer  than 
the  vacillation  period  and  from  six  to  nine  times  longer  than  the 
theirmal  lag  times  of  the  fluid.  The  shortest  imposed  temperature 
contrast  period  (experiment  SE-31)  was  7.8  times  as  long  as  the 
observed  vacillation  period  in  that  experiment  and  3 times  longer 
than  the  thermal  lag  time. 

(c)  The  thermal  response  times ; The  response  times  of  the 
zonally  averaged  temperature  (T)  at  different  radii  and  depths  to 
a step-change  in  external  conditions  depend  upon  the  physical 
properties  of  the  system  (e.g.,  mass,  thermal  diffusivity) , the 
efficiency  of  the  eddy  and  symmetric  convection  and  the  position 
relative  to  the  boundaries.  In  our  experiments  the  response 
times  to  a step-change  in  T.  were  observed  to  increase,  and  the 
magnitude  of  the  changes  in  T to  decrease,  with  distance  from 
the  cold  inner  wall.  For  exeonple,  in  experiment  SE-32  (series  F) 
the  e-folding  response  times  at  mid-depth  at  1/6,  1/2,  and  5/6 


n 
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gap-widths  measured  from  the  inner  wall  were  approximately  1140, 
1310,  and  1600  sec  (312,  358,  and  437  rot) , respectively.  The 
changes  in  T at  these  radii  were  approximately  1.8,  1.6,  and 
1.4®C,  respectively.  Moreover,  since  the  'step'  change  in  T* 
was  actually  accomplished  at  the  thermal  controls  of  the  circu- 
lating bath,  also  exhibited  a response  time.  In  experiment 
SE-32  this  time  was  measured  to  be  315  sec  (86  rot) , which  should 
be  subtracted  from  the  response  times  in  the  fluid  interior  in 
order  to  measure  the  effect  of  the  fluid  convection. 

(d)  The  lag  times:  The  mean  lag  times  and  amplitudes  of  the 
T variations  at  different  locations  in  the  annulus  in  response 
to  cyclic  variations  of  T^  depend  upon  the  same  factors  as  do  the 
thermal  response  times  and  amplitudes  following  step-changes  in 
T,.  In  addition,  they  depend  upon  the  frequency  of  the  imposed 
cyclic  variation.  In  our  experiments  the  mean  lag  times  were 
observed  to  increase,  and  the  amplitudes  of  the  T variations  to 
decrease,  with  distance  from  the  inner  wall.  With  progressively 
longer  periods  of  the  imposed  cyclic  variation  the  cunplitudes  of 
the  internal  temperature  variations  were  progressively  larger, 
and  the  lag  times,  although  longer  in  absolute  time,  represented 
a progressively  smaller  fraction  of  the  imposed  cycle.  These 
properties  may  be  illustrated  in  the  case  of  experiments  SE-05 
and  SE-29  (series -A)  in  which  T - was  cycled  with  periods  of 
984.1  and  1968  rot,  respectively,  between  temperature  extremes 
which  were  not  too  different  from  each  other.  The  amplitudes  of 
the  T variations  at  mid-depth  at  various  radii  in  the  annulus 


14 


were,  on  average,  smaller  by  50%  and  20%,  respectively,  than  the 
difference  in  amplitude  between  the  temperatures  at  the  same 
locations  in  null  experiments  corresponding  to  the  two  extremes 
of  the  imposed  cycle.  The  lag  times  in  these  experiments  are 
given  in  Table  3. 

6.  Low  pass  filtering  of  the  data 

Our  interest  in  the  present  study  is  primarily  in  the 
responses  of  the  fluid  system  to  imposed  variations  of  T,- . As 
indicated  in  Table  2,  the  periods  of  the  imposed  temperature 
variation  are  considerably  greater  than  the  periods  of  amplitude 
vacillation.  In  order  to  bring  out  more  clearly  the  longer- 
period  responses  to  the  imposed  temperature  variation,  we  used  a 
low-pass  filter  to  suppress  the  relatively  shorter-period  ampli- 
tude vacillation.  The  filter  was  constructed  in  the  following 
manner:  The  original  data  obtained  at  any  one  location  consists 
of  one  temperature  measurement  every  four,  five,  or  six  rot 
depending  upon  the  rotation  rate  of  the  experiment.  A moving 
average  of  21  consecutive  temperature  measurements  was  performed 
first,  with  each  point  given  equal  weight.  The  resulting  aver- 
ages were  then  subjected  to  a second  21-point  moving  average. 

These  averages  were,  in  turn,  subjected  to  a third  21-point  moving 
average  to  complete  the  construction  of  the  filter.  Each  of  the 
smoothed  temperature  values  in  the  resulting  time  sequence  con- 
tains contributions  from  temperatures  measured  at  61  different 
times  weighted  as  shown  in  Fig.  4.  This  figure  also  shows  the 


TABLE  3 


Lag  tines  at  mid-depth  at  three  radii  in  percent  of  cycle  period 
Experiment  Radius  (gap  width) 

iZi 1/2 5/6 

SE-05  16% 

SE-29  10% 


18% 

20% 

11% 

12% 

nurAbcr  of  data.  po\nts 


Fig.  4 Distribution  of  relative  weights  used  in  smoothing  the  data.  Distributions 
ate  shown  for  three  successive  passes  of  the  21-point  moving  average, in  . 
which  each  point  is  given  equal  relative  weight. 
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effect  of  each  pass  of  the  21-point  moving  average  on  the  relative 
weights  assigned  to  each  temperature  in  the  time  sequence. 

The  effect  of  smoothing  the  data  of  experiment  SE-29 
(series  A)  is  shown  in  Figs.  5a,  and  5b.  Fig.  5a  shows  the  time 
variation  of  both  the  smoothed  and  unsmoothed  square  of  the 
amplitude  of  wave  number  4 for  two  externally  imposed  cycles. 

Fig.  5b  shows  a similar  comparison  for  the  eddy  thermal  variance 
over  the  same  two  cycles.  The  imposed  variation  of  T-  is  de- 
picted  by  the  dotted  curve.  In  this  experiment,  in  which  the 
imposed  cycle  period  was  2 hours  (1968  rot) , and  the  basic  rota- 
tion rate  was  1.718  sec  the  low-pass  filter  represents  a 
weighted  average  over  244  rot,  which  is  approximately  1/8  the 
period  of  the  imposed  cycle.  In  most  of  our  experiments  the 
imposed  cycle  period  was  1 hour  (984.1  rot)  and  the  rotation  rate 
was  1.718  sec  in  which  case  the  low-pass  filter  represents  a 
moving  weighted  average  over  1/4  the  period  of  the  imposed  cycle. 
This  might  be  thought  of  as  a moving  weighted  "seasonal"  average 
with  the  heaviest  weighting  assigned  to  the  middle  "month" . 

7.  Interior  radial  and  vertical  temperature  differences 

The  low-pass  filter  was  applied  also  to  the  data  obtained 
from  experiments  in  which  the  thermostat  controlling  T^  was 
changed  in  steps.  The  responses  of  the  filtered  T's  at  mid-depth, 
at  1/6,  1/2,  and  5/6  gap-widths  measured  from  the  inner  wall, 
following  a near-step-change  in  T.  from  24.4  to  19.7®C  in  experi- 

V 

ment  SE-32  (series  F) , are  shown  in  Fig.  6. 
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Fig.  6;  Time  responses  of  the  inner  colc’-bath  temperature  and  of  the  smoothed  zonally- 
averaged  temperatures  at  mid-depth,  at  1/6,  1/2  and  5/6  gap-width  from  the 
inner  cylinder,  due  to  a step  change  of  the  thermoregulator  controlling  the 
inner-bath  temperature  in  experiment  SE-32. 


time  C 10^  ) 


24 

3 

23  " 


£ 

.3 


22 


21 


O 

X) 

u 

c 


20 


Fig.  7:  Time  responses  of  the  smoothed  zonal ly-averaged  temperature  at  mid-depth,  at 
1/6,  1/2  and  5/6  gap-width  from  the  inner  cylinder,  due  to  the  cyclic  varia- 
tion of  the  cold  inner-bath  temperature  in  experiment  SE-05 
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response  of  T-  to  the  change  in  the  thermostat.  The  scale  for 

Z' 

T • appears  on  the  right  hand  side  of  the  diagram.  The  inner-bath 
temperature  and  the  filtered  temperatures  in  the  interior  of  the 
working  fluid  appear  to  have  exponential  responses.  These  are 
the  curves  which  provided  the  data  for  the  e-folding  times  re- 
ported earlier. 

Fig.  7 shows  the  response  of  the  filtered  T ' s at  mid-depth, 
at  1/6,  1/2,  and  5/6  gap-widths,  to  cyclic  variations  of  T^ 
between  20.0  and  24.7®C  in  experiment  SE-05  (series  A).  The 
scale  for  T^  appears  on  the  right.  These  curves  provided  the 
data  for  the  lag  times  reported  earlier. 

From  the  data  in  Figs.  6 and  7 we  may  calculate  the  response 
of  the  filtered  temperature  difference  which 

serves  as  a crude  measure  of  the  zonal ly-averaged  temperature 
gradient  in  the  vicinity  of  mid-radius.  Since  the  zonally- 
averaged  temperatures  at  5/6  and  1/6  gap-width  have  different 
response  times  to  a step  change  in  T^ , due  to  their  different 
positions  relative  to  the  forcing  boundary,  the  response  curve 
for  the  difference  between  these  two  temperatures  is  not  a 
simple  decaying  exponential.  To  Illustrate  this  point  we  may 
consider  an  idealized  model  of  exponentially  decaying  temperature 
responses  at  the  two  radii 

Vs  = “ ^ 


d + ge 


-t/h 
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whereby  it  follows  that 


^5/6  - “^l/e  “ f + ge 


-t/h  _ be"^/^ 


(f  H d - a) 


(3) 


The  solid  curve  in  Fig.  8a  shows  the  idealized  response  curve 
for  the  difference  r using  the  values  a = 26.28®C, 

b = 1.80OC,  c = 1140  sec,  d = 27.22®C,  g = 1.44®C  and  h = 1600 
sec,  corresponding  to  experiment  SE-32  (Fig.  6).  The  most  notable 
feature  of  this  curve  is  the  considerable  overshoot  around  2000 
sec.  This  curve  may  be  compared  with  the  dashed  curve  which 
shows  the  response  which  would  occur  if  the  exponential  decay 
rate  of  T^^g  were  equal  to  that  of  if  h = c = 1140 

sec) . The  actual  response  of  the  difference  (T^^g  - 
experiment  SE-32,  shown  in  Fig.  8b,  displays  a steeper  and  larger 
overshoot  which,  as  shown  later  in  Fig.  13,  is  accompanied  by  a 
wave  number  change  at  the  time  of  the  peak.  Figure  9 shows  the 
observed  response  of  the  filtered  temperature  difference 
(Tg^g  - T^/g)  to  wave  number  transition  at  effectively  constant 
imposed  temperature  contrast  during  a section  of  experiment  SE-34 
(series  F) . The  notable  feature  here  is  the  peak  in  the  temper- 
ature difference  at  the  time  of  the  wave  number  change. 

We  consider  now  the  lag  time  of  the  smoothed  temperature 
difference  Tg^g  - ^2./6  ^ cyclic  variation  of  T^.  The 

temperature  variations  at  the  two  radii  may  be  expressed  in 
the  form 


t\me  C\o^sec) 


time  ( 10  rotations') 

Fig.  8a:  Theoretical  response  of  the  smoothed  tonally- averaged  temperature  difference 
(^5/6  “ yg)  to  a step-change  in  the  cold  bath  temperature  obtained  from 
equation  (3).  The  solid  curve  is  determined  using  coefficients  correspond- 
ing to  the  smoothed  mid-depth  data  of  experiment  SE-32.  The  dashed  curve 
is  obtained  using  the  same  values  of  f,  g,  b and  c,  but  setting  h=c=1140  sec 
corresponding  to  a hypothetical  case  in  which  the  exponential  decay  rate  of 
tsyg  is  equal  to  that  of  fiyg* 


-time , CIO*  sec.) 


time,  ( 10*  rotations) 


Observed  response  of  (ts.g  - Tl/s)  mid-depth  in  the  fluid  to  a near  step 
change  in  the  cold  inner-bath  tfcmperature  in  experiment  SE-32. 


18  19  20  Z\  Zt  Z^ 


time,  ( lo*  rotat\<jn<i) 


Observed  response  of  the  smoothed  temperature  difference  - Tj/g) 

associated  with  a change  in  the  flow  field  from  wave  number  4 to  wave 
number  5 at  effectively  a constant  imposed  ten^>erature  contrast  during 
experiment  SE-34. 
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Tj^^g  = m + n sin  (ot  - 

e) 

(4) 

Tgyg  = p + q sin  (ot  - 

n) 

(5) 

where  o is  the  frequency  of  the  oscillation  and  e and  n are  the 
time  lags  between  the  cold  bath  temperature  variation  and  the 
variations  of  ^5/6'  ^^®spectively.  It  follows  that 


Tgyg  - = r + q sin  (at  - n)  - n sin  (ot  - e)  , 


(r  s p - m) 


(6) 


or, 


^5/6  ■ Ve  = r - o) 


(7) 


where  D 


- * V * n^)  - 


) - 2nq  cos  (n  “ e)  , a 


- tan~^  { ^ ^ sine  \ 

\ q COST)  - n cose  / 


and  the  sign  of  D is  chosen  such  as  to  make  the  preceeding  two 
expressions  mutually  consistent.  The  solid  curve  in  Pig.  10a 
shows  the  idealized  response  curve  for  the  difference 
using  the  values  m = 26.71®C,  n = 0.47®C,  p 27.56®C,  q = 0.33®C, 
e = 56  deg  and  n = 71  deg  corresponding  to  experiment  SE-05 
(Fig.  7).  The  dashed  curve  shows  the  response  which  would  result 
if  n®e=56  deg.  The  notable  feature  here  is  the  larger  eunplitude 
and  earlier  response  of  the  solid  curve  to  the  imposed  variation 
of  T • (dotted  curve) . The  actual  response  of  the  smoothed 


0 0.2  0.^  ae  0.B  i.o 

time  ( 10^  rotatvo»'s) 


Fig.  10a:  Theoretical  response  of  the  smoothed  zonal ly-averaged  temperature  difference 
(ts.g  - Tl/e)  ^ cyclic  variation  of  the  cold  inner-bath  temperature  ob- 
tained from  equation  7.  The  solid  curve  was  determined  using  coefficients 
associated  with  the  smoothed  mid-depth  data  of  experiment  SE-05.  The  dashed 
curve  is  obtained  using  the  same  values  of  m,  n,  p,  q and  e,  but  setting 
n = e = 56®  corresponding  to  a hypothetical  case  in  which  the  lag  of  ts.e 
is  the  same  as  but  the  amplitude  of  is  smaller  than  that  or 

Tj  .g.  The  dash-dot-dashed  lines  represent  tne  estimated  maximum  and  minimum 
values  of  the  temperature  difference  (ts  .g  “ Ti /s)  corresponding  to  the 
maximum  and  minimum  imposed  temperature  contrast,  respectively,  in  the 
absence  of  cyclic  variation  of  the  cold  inner-bath  temperature. 


time,  ( 10*  sec.) 

0 lo  to 


Fig.  10b:  Observed_ response  of  the  smoothed  zonal ly-averaged  temperature  difference 

(T5/6  “ T) yg)  mid-depth  in  the  fluid  to  the  cyclic  variation  of  the  cold 
inner-bath' temperature  in  experiment  SE-05. 
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temperature  difference  'to  the  imposed  cyclic  varia- 

tion of  the  cold  bath  temperature  in  experiment  SE-05  is  shown  in 
Fig.  10b.  The  average  lag  between  the  minimum  T.  (i.e.,  maximum 
imposed  AT)  and  the  maximum  internal  temperature  difference 
(Tg^g  - Tj^yg)  is  30±9degvs  the  predicted  value  of  o = 26.5  deg 
shown  in  Fig.  13a. 

Finally,  we  present  the  response  of  the  smoothed  vertical 
temperature  difference  in  the  interior  of  the  fluid  to  the  imposed 
cyclic  variation  of  T^  in  experiment  SE-05.  The  network  of  ther- 
mistcrs  used  in  this  study  enabled  us  to  estimate  the  vertical 
stratification  only  at  1/6  gap-width  from  the  cold  wall  and  only 
in  the  layer  1/2  H < z < 3/4  H.  The  time  variation  of  the  tem- 
perature difference  between  these  two  levels  at  1/6  gap-width  is 
shown  by  the  solid  curve  in  Fig.  11.  The  stratifications  at 
larger  radii  can  be  expected  to  display  greater  time  lags  with 
respect  to  the  cyclic  variation  of  T . (dashed  curve) . 

V 

8.  Null  experiments 


Null  exper.  nts  are  defined  here  as  experiments  in  which 

the  externally  c .tcrolled  pareimeters  (rotation  rate  and  Imposed 

temperature  contrast)  are  held  fixed.  These  are  conventional 

annulus  experiments.  In  this  section  we  shall  discuss  pertinent 

null 

results  of  previously  unpublishedAexperiments  performed  with  a 
two-level  network  of  100  thermistors,  as  well  as  results  of  three 
additional  null  experiments  conducted  in  order  to  supplement  our 
knowledge  of  the  time-dependent  behavior  of  the  fluid  at 


time.  (10^  sec.) 


zz  y 


time  . ( 10^  rotcitloni) 


Observed  response  of  the  smoothed  vertical  temperature  diffemce  at  1/6 
gap-width  from  the  inner  cylinder  in  experiment  SE-05.  The  vertical 
temperature  difference  [(11,6)3.4  - (Ti,g)i/z1  was  obtained  from  the  zonally- 
averaged  temperatures  at  heights' z = 3/4' H and  z = 1/2  H. 
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significant  points  in  dimensionless-parameter  space.  Such  infor- 
mation provides  a useful  framework  within  which  to  interpret  the 
results  of  experiments  with  imposed  cyclic  variations  of  T • . 

(The  three  new  experiments  are  listed  in  Table  1 as  series  G and 
are  designated  in  Fig.  2 by  crosses.)  In  addition,  we  shall 
present  the  results  of  experiment  SE-32,  in  which  the  imposed 
AT  was  changed  abruptly  from  5.9  to  10.6®C  by  changing  T^  from 
24.4  to  19.7«C. 

In  all  null  experiments  presented  here  the  annulus  was 
maintained  stationary  at  first  while  the  hot  and  cold  tempera- 
ture baths  were  brought  to  the  required  equilibrium  values.  The 
rotation  rate  of  the  turntable  was  then  set  to  the  specified  j 

value.  After  initial  spin  up,  the  flow  normally  exhibited 
further  transient  characteristics  before  settling  down  to  a 
statistically  identifiable  behavior  which  we  define  as  the  char- 
acteristic behavior  at  the  particular  point  in  dimensionless- 
parameter  space. 

Figs.  12a,  b,  and  c present  contours  of  the  filtered  eddy 
temperature  variance  at  mid-depth  and  mid-gap-width  as  a func- 
tion of  wave  number  and  time  in  experiments  SE-37,  39,  and  25, 
respectively.  Contours  are  used  here  for  convenience  although 
it  is  recognized  that  the  variance  spectra  are  discrete.  Above 
each  figure  is  a plot  of  the  time  variations  of  the  filtered 
zonal  and  eddy  temperature  variances  (V^  and  V^,  respectively) 
evaluated  at  mid-depth,  and  of  the  imposed  cold  bath  temperatures. 

Experiments  SE-37,  39,  and  25  are  characterized  by  4-wave 
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amplitude  vacillation,  5-wave  amplitude  vacillation  and  4,  5,  3 
mixed  mode  amplitude  vacillation,  respectively.  In  the  4-  and  5- 
wave  amplitude  vacillation  cases  data  were  taken  starting  prior 
to  experimental  time  zero  (i.e.,  the  time  at  which  the  rotation 
rate  was  suddenly  changed  from  zero  to  the  specified  value) . 
Accordingly,  Figs.  12a  and  b show  clearly  the  transient  states 
in  these  two  experiments.  In  both  cases  the  transient  states 
are  charactertized  by  a broader  variance  spectrum  and  a larger 
zonal  temperature  variance  than  in  the  remaining  portion  of  each 
experiment.  It  is  interesting  to  note,  however,  that  although 
these  two  experiments  are  very  close  to  each  other  in  dimension- 
less-parameter space  (Ro^  = 0.81  and  Ro^  = 0.79,  respectively), 

( the  transient  is  shorter  in  duration  in  the  5-wave  amplitude 

vacillation  experiment  than  it  is  in  the  4-wave  vacillation 
experiment.  Another  interesting  feature  is  the  sharp  peak  in 
in  both  cases  just  prior  to  the  formation  of  the  long-term 
non- transient  pattern.  We  believe  this  is  due  to  the  ineffi- 
ciency of  the  broader  spectrum  of  waves  present  during  the 
transient  stage  in  transporting  heat  "poleward" , thus  leading 
to  a continued  rise  in  until  the  preferred  mode  becomes  well 
established  and  -begins  to  transport  sufficient  eunounts  of  heat 
"poleward"  to  balance  the  rate  of  generation  of  by  differen- 
tial heating.  As  noted  earlier,  amplitude  vacillations  are 
smoothed  out  of  the  data  by  the  low  pass  filter.  Accordingly, 


the  non-transient  portion  of  each  record  is  characterized  in 
Figs.  12a,  b,  and  c by  a steady  thermal  variance  maximum  at  the 
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wave  number  of  the  dominant  vacillating  mode,  with  much  smaller 
amounts  of  thermal  variance  at  higher  and  lower  wave  numbers. 

In  the  4,  5,  3 mixed  mode  amplitude  vacillation  experiment 
(Fig.  12c)  the  dominant  mode  is  clearly  wave  number  4,  but  there 
is  also  a definite  alternation  of  thermal  variance  with  time 
between  wave  numbers  3 and  4 on  the  one  hand  and  5 on  the  other. 
The  characteristic  time-scale  for  this  exchange  is  of  the  order 
of  125  to  150  rot.  This  is  several  times  longer  than  the 
amplitude  vacillation  period  of  wave  numbers  4 and  5,  but  about 
the  same  time  scale  as  the  fluctuations  of  wave  number  3 when  its 
amplitude  is  large. ■ff [Figure  13  presents  contours  of  the  filtered 
eddy  thermal  variance  as  a function  of  wave  number  and  time  in 
experiment  SE-32  in  which  the  imposed  temperature  contrast  was 
held  constant  at  AT  = 5.9®C  (the  same  value  as  in  SE-25,  Fig.  12c) 
for  more  than  8,000  rot  and  then  changed  precipitously  to  AT  = 
10.6®C  (which  is  not  too  far  from  the  value,  AT  = 9.0®C,  in 
SE-37,  Fig.  12a)  and  held  at  this  constant  value  for  more  than 
14,000  rot  before  being  returned  precipitously  to  AT  = 5.9®C. 
Above  the  figure  is  a plot  of  and  evaluated  at  mid-depth, 
and  of  T^.  At  each  value  of  the  imposed  AT  the  characteristic 
features  of  the  contour  patterns  and  energy  curves  are  similar 
to  those  displayed  in  null  experiments  at  the  same  or  nearby 
values  of  AT  (Figs.  12a  and  12c).  There  are  also  characteristic 
energy  levels  for  both  and  corresponding  to  each  value  of 
AT.  It  is  observed,  too,  that  when  T.  is  lowered,  around  t = 

X- 

8,150  rot,  resulting  in  a precipitous  increase  in  the  imposed  AT, 
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rises  sharply  to  a peak  emd  then  falls  to  a nev/  value  which  is 
larger  than  the  characteristic  value  at  the  lower  temperature 
contrast.  We  believe  that  the  sharp  increase  in  is  due  mainly  to 
the  inefficiency  of  the  "poleward"  eddy  heat  flux  associated 
with  the  pre-existing  flow  pattern  (in  this  case  mixed  mode 
amplitude  vacillation)  at  the  new  imposed  value  of  AT.  This 
enables  V^,  to  rise  in  response  to  the  precipitous  increase  in 
AT,  and  the  associated  internal  temperature  gradient  to  become 
baroclinically  unstable  for  the  formation  of  a new  flow  pattern 
(in  this  case  4-wave  amplitude  vacillation) . The  "poleward" 
eddy  flux  of  heat  associated  with  the  new  flow  pattern  is  be- 
lieved to  be  sufficiently  vigorous  to  reduce  Vj,  to  a new  char- 
acteristic level  consistent  with  this  flow  pattern.  It  will  be 
noted,  however,  that  the  reverse  change  in  AT  from  10.6  to 
S.S^C  (around  t = 22,570  sec)  is  not  accompanied  by  a similar 
j overshoot. 

I Results  derived  from  previously  unpublished  null  experiments 

' which  have  an  important  bearing  on  the  results  to  be  expected 

, from  experiments  with  cyclic  variations  of  T . were  presented 

earlier  in  Pig.  3.  In  particular,  this  figure  shows  the  relative 
magnitude  of  the  eddy  thermal  variance  as  a function  of  the  rota- 
tion rate  at  three  different  externally  imposed  temperature 
contrasts.  Here  we  see  that  the  rate  of  increase  of  eddy  thermal 
variance  with  rotation  rate  is  greater  at  successively  greater 
values  of  the  imposed  AT.  From  the  arrangement  of  the  numbers, 
we  can  anticipate  that  in  experimental  series  A and  E there 

i 


should  be  little  if  any  change  in  the  eddy  thermal  variance  from 
"summer"  to  "winter"  as  T.  is  varied  cyclically.  Experiments  C 
and  D,  on  the  other  hand,  would  be  expected  to  display  greater 
eddy  thermal  variance  in  "winter"  (i.e.,  at  the  higher  values  of 
the  imposed  temperature  contrast)  than  in  "summer",  as  does  the 
earth ' s atmosphere . 

Another  feature  of  interest  seen  in  Figs.  2 and  3 is  that 
two  different  wave  regimes  (i.e.,  5-wave  amplitude  vacillation 
and  6-wave  amplitude  vacillation)  occurred  when  the  experiment 
at  Ro^  =0.50  and  Ta  = 1.43  x 10^  was  run  on  different  days.  The 
higher  wave-number  flow  is  seen  to  possess  a lower  level  of  eddy 
thermal  variance.  A similar  result  occurs  in  other  experimental 
series  at  various  points  in  dimensionless-parameter  space. 

9.  Hysteresis  experiments 

Hysteresis  experiments  are  ones  in  which  one  of  the  exter- 
nally imposed  parameters  (e.g.,  AT  or  0)  is  varied  extremely 
slowly  from  an  initial  to  a final  value  and  then  back  to  the 
initial  value.  If  hysteresis  is  present  it  manifests  itself  as 
a bi-modal  response  of  the  system  to  the  forcing  function  such 
that  different  fluid  behavior  is  observed  at  a given  point  in 
dimensionless-parameter  space  when  it  is  traversed  from  different 
directions.  Hysteresis  in  rotating  annulus  convection  has  been 
studied  experimentally  by  Fultz,  et  al  (1959)  and  Kaiser  (1970), 
and  theoretically  by  Cole  (1971).  This  phenomenon  has  also  been 
studied  in  terms  of  catastrophe  theory  by  Lacher,  McArthur  and 
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Buzyna  (1977).  In  this  section  we  present  some  details  of  a 
hysteresis  study  conducted  in  the  region  of  the  regime  diagram 
in  which  experimental  series  A was  performed.  The  results  pro- 
vide background  information  which  is  useful  for  interpreting 
the  data  from  experiments  with  cyclic  variations  of  T . . 

Three  hysteresis  experiments  were  conducted  (series  F, 
experiments  SE-34,  35,  and  36).  In  experiment  SE-34  the  exter- 
nal temperature  contrast  was  varied  over  a range  of  4.5^C  by 
varying  the  cold  bath  temperature  in  four  steps  over  a period 
of  16  hours.  Since  the  cold  bath  temperature  at  each  step  was 
maintained  constant  for  2 hours,  the  thermal  response  time  of 
the  fluid  (which  is  about  1/5  this  length  of  time)  is  not  an 
important  factor  for  present  considerations.  Fig.  14  shows  con- 
tours of  the  smoothed  thermal  variance  as  a fvmction  of  wave 
number  and  time  in  this  experiment.  The  step-like  dashed  curve 
above  the  contour  diagram  represents  the  time  variation  of  T>. 

The  imposed  temperature  contrast,  AT,  indicated  at  each  step, 
was  changed  from  10.6®C  at  the  beginning  of  the  experiment  to 
6.1®C,  and  then  back  to  10.6®C.  The  experiment  was  continued  for 
another  half  cycle  to  AT  = 5. 8®C.  The  upper  and  lower  solid 
curves  above  the  contour  diagreun  represent  the  time  variations 
of  the  smoothed  and  V^,  respectively,  measured  at  mid-depth. 

Fig.  14  reveals  that  hysteresis  is  indeed  an  important 
characteristic  of  the  fluid  behavior.  The  experiment  begins 
with  AT  = 10.6®C,  at  which  value  wave  number  4 is  the  dominant 
mode.  (The  amplitude  vacillations  of  this  and  other  modes 


represented  and  contour  values  and  intervals  are  the  same  as  in  Figs  12  and  13.  Numbers  above  each  step  of 
the  cold  bath  temperature  indicate  the  values  of  the  imposed  temperature  contrast  corresponding  to  that  step 
The  vertical  stratification  used  in  calculating  the  energy  equivalent  of  V_  and  is  determined  from  the 
time  average  of  [(ti.6)3.4  - (tj, 6)1.2]  over  each  of  the  steps  in  the  first  complete  cycle. 


throughout  the  experiment  are  not  evident  in  the  diagram  because 
of  the  smoothing  function  applied  to  the  data.)  Later,  when  the 
imposed  temperature  difference  is  returned  to  the  value  AT  = 10.6 
°C  (15,000  < t < 16,890  rot)  a preexisting  5-wave  pattern  changes 
to  a predominantly  4-wave  pattern.  Apparently,  a 10.6®C  imposed 
temperature  contrast  can  sustain  a 5-wave  pattern  for  only  a 
limited  period  of  time,  after  which  a change  takes  place  to  the 
preferred  4-wave  pattern.  It  is  noteworthy  that  V2  peaks  in 
this  interval,  at  which  time  wave  number  5 ceases  to  be  the 
dominant  mode.  The  peak  around  t = 15,720  rot  is,  in  fact,  the 
largest  magnitude  attained  by  during  the  entire  experiment. 

At  this  point  it  appears  that  the  flow  pattern  becomes  baro- 
clinically  unstable  for  the  formation  of  a broader  spectrum  of 
wave  scales  of  which  wave  nunber  4 is  the  fastest  growing  mode. 
The  peak  in  apparently  denotes  the  time  at  which  the  vigor- 
ous "poleward"  heat  flux  accomplished  primarily  by  the  newly 

formed  4-wave  pattern  must  first  exceed  the  rate  of  generation 
. 3 

of  Vg  by  diabatic  heating. 

The  second  step  corresponds  to  an  imposed  temperature  con- 
trast of  9.8®C.  When  wave  number  4 is  present  at  the  beginning 

^In  speaking  of  the  "poleward"  heat  flux  accomplished  by 
one  or  more  wave  scales  we  must  recognize  that,  by  virture  of 
their  heat  and  momentum  fluxes,  different  waves  generate  differ- 
ent distributions  of  zonal ly- averaged  meridional  circulations 
which  oppose  the  heat  flux  in  certain  latitudes  and  enhance  it 
in  others.  The  efficiency  of  any  wave  pattern  in  accomplishing 
a "poleward"  heat  flux  must  be  assessed  in  terms  of  the  net 
effect  of  both  the  wave  pattern  itself  and  that  portion  of  the 
symmetric  meridional  circulation  which  it  generates. 
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it  remains  during  this  step.  However,  when  wave  number  5 is 
already  present,  as  it  is  when  the  condition  aT  = 9.8°C  is 
traversed  from  the  opposite  direction  (during  the  time  interval 
13,090  i t < 14,870  rot),  the  5-wave  configuration  remains 
during  the  entire  step.  Thus,  a 9.8°C  imposed  temperature  contrast 
is  capable  of  sustaining  either  a 4-  or  5-wave  amplitude  vacil- 
lating pattern  depending  upon  the  initial  conditions. 

The  third  step  corresponds  to  AT  = 8.2®C.  At  this  temper- 
ature contrast  the  preexisting  4-wave  cimplitude  vacillation 

comparative ly 

pattern  can  be  sustained  only  for  a^i short  period  of  time  before 
it  changes  to  a predominantly  5-wave  amplitude  vacillation 
pattern.  It  is  significant  that  there  is  a peak  in  at  the 
time  the  wave  number  changes,  suggesting  that  the  newly-formed 
5-wave  pattern  at  this  point  in  dimensionless-parameter  space 
is  more  efficient  at  transporting  heat  "poleward"  than  is  the 
preexisting  4-wave  pattern  which  permitted  to  rise  to  this 
peak.  Later,  when  the  imposed  temperature  contrast  aT  = 8.1°C 
is  traversed  from  the  opposite  direction,  during  the  time  inter- 
val 11,320  ^ t < 12,850  rot,  the  preexisting  5-wave  pattern  is 
sustained  throughout  the  entire  interval. 

The  fourth  step  corresponds  to  AT  = 6.7®C.  Here  we  see  a 
tendency  for  spectrum  broadening  and  a slight  wavering  of  the 
preexisting  5-wave  pattern.  Later  when  AT  = 6.6®C  is  traversed 
from  the  opposite  direction,  a preexisting  mixed  mode  amplitude 
vacillation  changes  to  a wavering  5-wave  amplitude  vacillation. 
Finally,  at  AT  = 6.1°C  the  weakly  wavering  dominant  5-wave 
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configuration  changes  to  a predominantly  mixed  mode  vacillation 
which  seems  to  be  the  preferred  configuration  at  this  value  of 
the  imposed  temperature  contrast.  The  final  third  of  the  figure 
displays  essentially  the  saune  behavior  under  essentially  the 
same  imposed  external  conditions  as  the  first  third.  In  addition 
to  the  characteristics  discussed  thus  fax;  it  is  also  noteworthy 
that  wave  number  changes  are  preceded  and  accompanied  by  spectral 
broadening.  It  may  be  recalled  that  this  experiment  was  conducted 
in  a region  of  the  regime  diagram  (Fig.  3)  in  which  the  eddy 
thermal  variance  is  no  greater  for  AT  = 10°C  them  it  is  for 
AT  = 5°C.  Hence,  the  only  observed  changes  in  the  magnitude  of 
Vg  (lower  solid  curve  in  the  figure)  are  those  associated  with 
wave-number  changes.  In  this  connection,  it  is  significant  that 
the  magnitude  of  is  measurably  lower  when  wave  number  5 is 
the  dominant  vacillating  mode  than  when  the  flow  is  characterized 
by  either  a predominantly  4-wave  or  mixed-mode  amplitude  vacilla- 
tion. 

We  turn  now  to  Fig.  15  which  displays  the  results  of  exper- 
iment SE-35.  In  this  experiment  the  imposed  temperature  contrast 
was  varied  in  smaller  steps  over  the  limited  range  8.1  ^ AT  i 10.5°C 
with  fl  set  at  the  same  value  as  it  was  in  experiment  SE-34.  In 
SE-35  the  cold-bath  temperature  at  each  step  was  maintained 
constant  for  3 ho’ar.s.  This  experiment  began  with  AT  = 8.1®C,  at 
which  a dominant  5-wave  cunplitude  vacillation  developed  and  was 
sustained  as  long  as  the  temperature  contrast  was  maintained  at 
this  value.  This  behavior  is  consistant  with  that  seen  in  Fig.  14. 
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Later,  when  the  same  imposed  temperature  contrast  was  approached 
from  higher  values  of  AT,  a preexisting  mixed-mode  pattern 
changed  to  a 5-wave  amplitude  vacillation.  In  experiment  SE-34 
(Fig.  14),  when  the  temperature  contrast  AT  = 8.0°C  was  approach- 
ed from  AT  = 9.8°C  (representing  a larger  jump)  there  was  a 
change  from  a 4-wave  to  a 5-wave  eunplitude  vacillation  with  a 
brief  period  of  spectral  broadening  separating  the  two  wave 
states.  It  appears  from  excimination  of  these  two  experiments 
that  the  broader  spectral  energy  distribution  is  the  preferred 
condition  when  AT  is  being  decreased  slowly  in  the  interval 
8.1  < AT  < 8.5°C.  When  this  interval  is  traversed  in  small 
stepsr  the  mixed-mode  behavior  is  sustained  (Fig.  15),  but,  when 
it  is  traversed  by  a large  stepy  the  mixed-mode  behavior  appears 
only  briefly  between  4-  and  5-wave  amplitude  vacillation.  These 
results  suggest  that  the  fluid  behavior,  at  least  in  certain 
regions  of  dimensionless-parameter  space,  depends  upon  both  the 
direction  and  the  speed  of  approach  from  neighboring  regions. 

It  should  be  noted,  however,  that  certain  of  the  steps  in  Figs, 
either 

14  and  15  displayAunique  or  bi-modal  behavior.  For  exeimple, 
in  Fig.  14  we  see  that  when  AT  = 6.7°C,only  5-wave  amplitude 
vacillation  was  observed,  and  when  AT  ■»  9. 8®C,  either  4-wave  or 
5-wave  amplitude  vacillation  was  observed,  depending  upon 
whether  AT  was  being  decreased  or  increased. 

The  third  step  in  Fig.  15,  corresponding  to  AT  *»  9.2®C 
produces  a slight  change  in  the  5-wave  amplitude  vacillation 
pattern.  In  particular,  the  eunplitude  of  wave  number  5 decreases 
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and  the  spectrum  tends  to  broaden  slightly  toward  lower  wave 
numbers,  suggesting  that  S-wave  eunplitude  vacillation  is  not  as 
well  established  as  it  was  in  the  preceding  step.  When  a temper- 
ature contrast  close  to  this  value  (AT  = 9.3°C)  is  traversed  from 
the  opposite  direction,  a preexisting  4 -wave  eimplitude  vacilla- 
tion is  sustained  during  the  entire  step. 

The  fourth  step,  corresponding  to  AT  = 10.1°C,  produces  a 
change  from  the  5-wave  to  a 4-wave  eunplitude  vacillation  in  the 
interval  8,260  i t i 11,130  rot;  but  much  later  during  the 
second  cycle  (32,130  i t ^ 34,830  rot)  the  preexisting  5-wave 
amplitude  vacillation  is  only  weakened  and  the  spectrum  broaden- 
ed slightly  toward  lower  wave  numbers  when  AT  is  changed  from 
9.1  to  10.0°C.  When  a temperature  contrast  close  to  this  value 
(AT  s=  10.2°C)  is  traversed  from  the  opposite  direction  during 
the  time  interval  14,220  1 t 1 17,090  rot  a preexisting  4-wave 
amplitude  vacillation  is  seen  to  persist.  The  final  step, 
corresponding  to  AT  = 10.5®C  sustains  the  preexisting  4-wave 
amplitude  vacillation  during  the  first  cycle  and  changes  the 
5-wave  to  a 4-wave  amplitude  vacillation  during  the  second  cycle. 
These  results  are  consistent  with  those  of  Fig.  14. 

In  addition  to  the  two  experiments  described  above,  we  also 
ran  an  experiment  (SE-36)  in  which  AT  was  varied  slowly  and 
continuously  over  essentially  the  same  range  as  in  SE-35  with  a 
24-hour  cycle  period.  Unfortunately,  difficulties  with  our  data 
system  created  large  data  gaps  in  the  record.  Nevertheless,  it 
was  considered  worthwhile  to  present  the  results  (Fig.  16) 
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because  of  their  consistency  with  those  shown  in  Fig.  15. 

Fig.  17  is  a schematic  diagram  encompassing  the  results  of 
the  hysteresis  experiments  presented  in  Figs  14,  15,  and  16.  We 
may  summarize  these  results  as  follows:  When  (2  is  held  constamt 
at  1.718  sec  ^ and  AT  is  changed  sufficiently  slowly  that  the 
fluid  has  ample  time  to  adjust  to  each  new  imposed  temperature 
contrast,  5-wave  amplitude  vacillation  is  the  dominant  character- 
istic in  the  interval  6.6  < AT  < 8.0°C  regardless  of  whether  AT 
is  being  increased  or  decreased  through  this  interval.  The  5- 
wave  aimplitude  vacillation  will  sustain  itself  if  AT  is  increased 
in  the  interval  8.0  < AT  < 10.1®C,  but  will  change  to  a dominant 
4-wave  amplitude  vacillation  when  AT  exceeds  10.1°C>.  The  4-wave 
amplitude  vacillation  will  prevail  if  AT  is  brought  back  in  steps 
to  8.5®C  and  will  change  to  a mixed  mode  amplitude  vacillation  in 
the  interval  8.5  > AT  > 8.0°C  before  changing  to  the  dominant  5- 
wave  6unplitude  vacillation  which  is  characteristic  of  the  inter- 
val 6.6  < AT  < 8.0°C.  When  AT  is  lowered  further  the  5-wave 
amplitude  vacillation  persists  until  AT  6.1°C  at  which  time  a 
mixed  mode  aunplitude  vacillation  develops.  As  AT  is  raised 
again  this  pattern  is  preserved  until  AT  6.6°C  at  which  time 
the  pattern  returns  to  a dominant  5-wave  zunplitude  vacillation. 
Thus,  the  history  of  the  flow  determines  the  observed  flow 
pattern  in  the  range  6.1  < AT  <10.6°C,  except  within  the  interval 
6.5  < AT  < 8.0®C  where  wave  number  5 is  observed  to  dominate 
whether  AT  is  being  increased  or  decreased.  Similar  behavior 
can  be  expected  in  other  regions  of  dimensionless-parameter 
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Fig.  17:  Schematic  diagram  summarizing  the  results  of  the  hysteresis  experiments 

SE-34,  35  and  36  shown  in  the  preceding  three  figures.  The  dominant  modes 
present  in  the  flow  field  are  indicated  as  a function  of  the  imposed 
temperature  contrast,  AT,  and  its  direction  of  change. 
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space. 

Fig.  17  displays  three  regions  of  mixed  mode  amplitude 
vacillation  which  appear  to  separate  different  regions  of  pure- 
mode amplitude  vacillation.  These  are  6.1  < AT  < 6.6°C, 

8.0  < AT  < 8.5°C  and  10.1  < AT  < 10.6°C.  Close  examination  of 
the  data  corresponding  to  each  of  these  regions  reveals  that 
the  characteristic  behavior  of  the  fluid  is  different  in  the 
interval  6.1  < AT  < 6.6°C  than  it  is  in  the  other  two  intervals. 

The  characteristics  of  the  contour  diagram  depicting  the  thermal 
variance  as  a function  of  wave  number  and  time  in  this  region 
are  similar  to  those  in  Figs.  12c  and  13,  where  AT  = 5.9°C.  In 
particular,  there  appears  to  be  an  alternation  •'f  thermal  vari- 
ance between  wave  number  5 on  the  one  hand  and  wave  numbers  3 and 
4 on  the  other.  (As  noted  earlier  in  connection  with  Fig.  12c, 
the  time-scale  for  this  exchange  is  several  times  longer  than 
the  amplitude  vacillation  period)  . On  the  other  hand,  the  spec- 
tral broadening  in  the  other  two  intervals  which  takes  place 
either  briefly  or  for  a sustained  period  of  time  during  wave- 
number  changes  from  4 to  5 or  from  5 to  4 displays  the  simul- 
taneous presence  of  all  three  modes. 

10.  Cyclic  variations  of  the  imposed  temperature  contrast 

Ten  experiments  were  conducted  with  cyclic  variations  of 
the  Imposed  temperature  contrast.  These  are  included  within  series 
A,  B,  C,  D,  and  E in  Table  1 and  in  Fig.  1.  In  each  experiment 
a null  state  was  established  first,  corresponding  either  to  the 
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maximum  or  minimum  temperature  contrast  in  the  cycle.  Next,  the 
temperature  contrast  was  varied  cyclically  between  specified 
limits.  The  shortest  experiment  comprised  5 cycles;  the  longest, 
50  cycles.  Cycle  periods  used  were  2,  1 and  0.5  hours,  covering 
a range  of  periods  from  2516  to  492.1  rot.  The  lengths  of  the 
cycles  were  from  8 to  57  times  longer  than  the  amplitude  vacil- 
lation periods  (see  Table  2)  and  from  4 1/2  to  9 times  longer 
thermal 

than  the  average/flag  time  of  the  fluid  at  mid-radius. 

In  interpreting  the  results  of  these  experiments  we  shall 
make  frequent  reference  to  Fig.  17.  This  figure  reveals  that, 
when  AT  is  changed  slowly  in  the  range  6.1  5 AT  i 10.6  C,  the 
possible  modes  of  behavior  in  each  of  five  intervals  are  limited 
and  predictable  as  a function  of  AT.  With  the  use  of  defini- 
tions introduced  into  the  meteorological  literature  by  Lorenz 
(1968)  , the  behavior  of  the  fluid  would  be  characterized  as 
transitive  within  the  interval  6.6  < AT  < 8.0°C  and  intransi- 
tive outside  this  interval.  A flow  is  said  to  be  transitive  if 
its  long-time-averaged  statistics  are  independent  of  the  initial 
state,  and  intransitive  if  its  long-time-:averaged  statistics 
depend  upon  the  initial  state.  In  the  interval  6.6  < AT  < 8.0®C, 
five-wave  amplitude  vacillation  occurred  whether  the  imposed 
temperature  contrast  was  being  increased  or  decreased.  Outside 
this  interval,  the  observed  flow  configuration  depended  upon  the 
direction  from  which  the  imposed  temperature  contrast  was 
approached. 

Both  transitive  and  intransitive  flows  are  determinate  in 
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the  sense  that  both  are  predictable  once  the  finite  amplitude 
flow  becomes  sufficiently  well  established.  When  AT  is  cycled  more 
rapidly  through  the  same  range,  however,  the  fluid 

may  not  have  adequate  time  to  respond  fully  to  the  changing 
conditions. 


A well-established  flow  configuration  might,  therefore,-^ 

be  sustained  through  portions  of  the  cycle  in  which  the 
imposed  conditions  are  otherwise  unfavorable  for  its  development. 
Moreover,  a flow  pattern  which  is  poorly  established  because  of 
the  rapidity  with  which  the  imposed  conditions  passed  through  the 
interval  most  favorable  for  its  development  might  easily  be  upset 
and  replaced  by  another  pattern  even  before  reaching  the  limit 
indicated  in  Fig.  17  for  the  onset  of  that  pattern.  As  a result 
of  such  possibilities,  the  fluid  could  display  different  behav- 
ioral characteristics  at  the  same  stage  of  different  imposed 
cycles,  depending  upon  how  well  established  the  pattern  is  as  it 
enters  each  stage.  We  shall  characterize  such  behavior  as 
indeterminate . We  shall  define  a flow  as  indeterminate  if,  by 
virtue  of  its  own  internal  non-linear  dynamics,  -he  time  series 
of  any  significant  parameter  describing  the  flow’ can  alternate 
sporadically  between  two  or  more  characteristic  flow  configura- 
tions. "Almost  intransitive"  flow,  as  defined  by  Lorenz  (1968), 
can  be  considered  to  be  a special  case  of  indeterminate  flow  in 
which  given  flow  configurations  last  for  long  time  periods  in 
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comparison  with  the  characteristic  time  scales  of  the  fluid. 

On  average,  it  is  reasonable  to  expect  that  the  faster 
we  cycle  the  imposed  temperature  difference  through  the  pre- 
scribed range,  the  greater  will  be  the  lag  between  it  and  the 
internal  gradient  of  the  zonally-averaged  temperature  (measured 
as  a percentage  of  the  cycle  time) , and  the  smaller  will  be  the 
excursion  of  the  internal  temperature  gradient,  3T/3R,  from  its 
time-mean  value.  In  a particular  phase  of  a particular  cycle, 
however,  the  actual  radial  temperature  distribution  in  the  inter- 
ior of  the  fluid  will  depend  also  upon  the  nature  of  the  convec- 
tive processes  present  during  that  phase  of  the  cycle  under 
consideration.  For  example,  when  AT  is  intensified,  3T/3R  will 
tend  to  increase  faster  in  the  presence  of  a wave  pattern  which 
is  inefficient  in  transporting  heat  "poleward"  than  it  would  in 
the  presence  of  one  which  is  more  efficient.  In  an  indeterminate 
system,  3T/3R  could  be  strong  or  weak  at  the  same  phase  of  different 
imposed  cycles, depending  upon  which  wave  scale  dominates  and 
how  well-established  it  is  when  it  enters  this  phase  of  each  cycle. 
In  the  presence  of  an  inefficient  wave  patteriv  3T/3R  could  in- 
crease to  the  point  of  baroclinic  instability  for  the  formation 
of  a new,  more  efficient,  wave  pattern,  after  which  it  would 
decrease  more  rapidly  than  would  be  expected  from  simple  argu- 
ments concerning  the  response  of  the  fluid  to  the  cyclic  vari- 
ation of  AT.  Irregularities  in  the  behavior  of  3t/3r  in  the 
interior  of  the  fluid  that  would  be  implied  by  such  behavior 
were  observed  in  the  data  (see,  for  example.  Fig.  10b). 
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Because  of  the  strong  relationship  between  3T/3R  and  V„ , the 

z 

existence  of  such  irregularities  is  reflected  also  in  the  curves 
depicting  the  time  variation  of  to  be  presented  in  this  sec- 
tion. 

We  consider  first  the  experiments  in  series  A,  all  of  which 
were  conducted  at  = 1.718  sec”^.  The  experiment  with  the 
longest  cycle-period  in  this  series  (viz.,  1968  rot)  was  SE-29. 
Fig.  18  shows  contours  of  thermal  variance  as  a function  of 
wave  nximber  and  time,  and  curves  of  and  V^.  as  a function  of 
time,  corresponding  to  this  experiment.  Initially,  the  cold- 
bath  temperature  was  sustained  at  its  highest  value  in  this  ex- 
periment, producing  the  minimum  imposed  temperature  contrast, 

5.9®C.  The  flow  during  this  time  interval  was  an  alternating 
mixed-mode  amplitude  vacillation  similar  to  that  found  in  SE-25 
(Fig.  12c)  and  in  the  first  part  of  SE-32  (Fig.  13) , during 
which  AT  = 5.9®C.  Such  behavior  was  also  observed  in  the  hystere- 
sis experiment,  SE-34  (Fig.  14),  while  the  imposed  temperature 
contrast  was  undergoing  intensification  in  the  interval 
6.1  5 AT  < 6.6«C. 

After  the  initial  sustained  "summer"  condition,  the  imposed 
temperature  difference  between  the  cold  and  hot  baths  was 
cycled  in  the  range  5.9  1 AT  1 10.6®C.  The  contour  diagram  for 
this  experiment  displays  at  different  times  the  various  char- 
acteristics of  the  temperature  spectra  observed  in  the  hystere- 
sis experiments  and  summarized  in  Fig.  17.  In  particular,  the 
diagreun  reveals  time  intervals  with  predominantly  5-wave  and 
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4-wave  aunplitude  vacillation  broken  by  briefer  intervals  of  mixed- 
mode amplitude  vacillation.  Unlike  the  hysteresis  experiments, 
however,  corresponding  phases  of  the  imposed  cycle  do  not  neces- 
sarily yield  the  same  temperature  spectriim.  For  exeunple,  5-wave 
amplitude  vacillation  sometimes  persists  through  one  or  two  cycles, 
and  4-wave  amplitude  vacillation  appears  less  frequently  and  for 
shorter  durations.  Although  every  "winter"  is  not  accompanied 
by  4-wave  £unplitude  vacillation,  every  appearance  of  4-wave 
amplitude  vacillation,  however,  occurs  in  "winter"  (i.e.,  in  the 
quarter  cycle  beginning  with  the  lowest  value  of  the  cold  bath 
temperature).  Similarly,  although  every  "summer"  is  not  accom- 
panied by  alternating  mixed-mode  amplitude  vacillation,  prac- 
tically every  appearance  of  this  form  of  vacillation  occurs  in 
"summer".  As  we  might  have  anticipated  from  Pig.  17,  5-wave 
cimplitude  vacillation  is  capable  of  persisting  through  the 
extremes  of  "winter"  and  "summer"  when  this  flow  configuration 
is  well  established,  and  when  the  ratio  of  the  imposed  cycle 
period  to  the  thermal  lag  time  is  much  smaller  than  it  is  in  the 
hysteresis  experiments.  In  experiment  SE-29  this  ratio  is  8.4 
which  is  of  the  order  of  10  times  smaller  than  it  is  in  the 
hysteresis  experiments.  Experiment  SE-29  clearly  displays 
indeterminate  behavior.  It  will  be  noted,  too,  that  there  are 
irregularities  in  the  time  variation  of  whenever  wave-number 
changes  take  place  and  that  the  magnitude  of  V^.  is  consistently 
smaller  when  5-wave  zunplitude  vacillation  dominates. 

It  might  be  expected  that  experiments  conducted  with  a 
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higher  frequency  of  cycling  within  the  same  limits  of  imposed  AT, 
and  at  the  same  value  of  fl,  would  display  even  more  erratic 
behavior,  with  the  internal  non-linear  dynamics  of  the  flow 
exerting  greater  control  over  events  and  the  particular  phase 
of  the  imposed  cycle  playing  a smaller  role.  In  this  connection, 
we  may  examine  the  results  of  experiments  SE-23,  24,  and  26 
(Figs.  19,  20,  and  21,  respectively)  in  which  the  imposed  cycle 
period  was  984.1  rot,  with  Q = 1.718  sec  and  AT  was  varied 
over  the  intervals  6.1  1 AT  i 10.6°C,  5.9  i AT  i 10.4°C  and 
5.6  i AT  i 10.2°C,  respectively.  The  control  parcuneters  were 
sufficiently  close  to  justify  considering  these  experiments 
together  as  a group.  Although  our  data  acquisition  system  mal- 
functioned on  several  occasions  during  each  experiment,  leaving 
periods  of  missing  data,  there  is  an  ample  amount  of  information 
on  which  to  derive  conclusions.  The  most  striking  feature  ob- 
served upon  comparing  the  results  of  the  three  experiments  is 
the  fact  that,  although  5-wave  cimplitude  vacillation  is  inter- 
rupted by  other  modes  of  behavior  on  several  occasions  during 
experiments  SE-23  and  SE-26,  this  phenomenon  continues  uninter- 
rupted through  19  of  the  20  "annual"  cycles  in  experiment  SE-24. 
We  can  only  speculate  that  the  5-wave  amplitude  vacillation  in 
SE-24  was  more  firmly  established  than  was  the  case  in  the  other 
two  experiments  and,  in  spite  of  numerous  tendencies  toward 
spectrxim  broadening,  a condition  which  usually  precedes  wave 
number  changes,  was  never  sufficiently  unstable  for  the  forma- 
tion of  another  wave  pattern.  In  any  case,  it  is  evident  that 
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TIME  I »*  HBTATiOM*) 

Response  of  the  temperature  variance  to  1-hour  cyclic  varia^pns  of  the  cold  inner-bath  temperature 
over  a 4.6*^  range, providing  imposed  temperature  contras tSAttom  5.6  to  10.2®C  in  experiment  SE-26, 
series  A (n  » 1.718  sec-1).  The  remaining  comments  are  the  same  as  in  Figs.  18  and  19. 
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the  internal  non-linear  dyneunics  of  the  flow  plays  a highly  signi- 
ficant role  when  the  imposed  temperature  difference  is  cycled  as 
rapidly  as  it  was  in  these  experiments,  in  which  the  imposed  cycle 
period  was  5.5  times  as  long  as  the  thermal  lag  time.  The  fluid 
behavior  under  these  conditions  is  clearly  indeterminate.  In 
such  a situation  it  would  appear  that  imperceptible  differences 
between  two  experiments  with  essentially  the  same  imposed  condi- 
tions, or  between  two  or  more  parts  of  the  same  experiment,  deter- 
mine whether  wave-pattern  changes  will  take  place,  and  at  what 
points  in  time  such  changes  will  occur. 

It  will  be  noted  that  in  all  three  experiments  the  most 
pronounced  irregularities  in  reflect  wave  number  changes, 
and  that  is  systematically  smaller  when  5-wave  amplitude 
vacillation  is  taking  place  than  it  is  at  other  times.  It  is 
also  noteworthy  that  when  the  cycling  is  terminated  with  an 
extended  "winter"  or  "summer"  condition  at  the  end  of  SE-23 
and  SE-24  the  temperature  spectrum  adopts  the  appropriate  con- 
figuration for  that  season. 

We  next  examine  Figs.  22  and  23  depicting  the  results  of 
experiments  SE-04  and  SE-05  in  which  the  imposed  cycle  time  and 
rotation  rate  were  again  984.1  rot  and  1.718  sec”^,  respectively. 

The  range  of  imposed  temperature  differences  in  both  experiments 

was  5.2  i AT  i 9.9°C,  the  mean  of  which  is  somewhat  lov.’er  than 

those  in  SE-23,  24,  and  26.  Experiment  SE-04  begins  with  a 

sustained  "winter"  condition  and  a 4-wave  funplitude  vacillation,  « 

but  changes  to  a flow  pattern  consisting  of  an  oscillating 
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mixed-mode  cunplitude  vacillation  characteristic  of  the  "summer" 

flows  and  a broad  spectrxim  characteristic  of  wave  form  changes. 

Experiment  SE-05  begins  with  a sustained  summer  condition  and  has 

characteristics  similar  to  those  in  SE-04.  Although  the  range  of 

imposed  temperature  differences  covers  the  region  of  5-wave  eunpli 

tude  vacillation  in  Fig.  17,  this  phenomenon  was  not  observed 

for  any  sustained  period  of  time  during  either  experiment.  In 

the  course  of  alternations  of  variance  between  wave  nximber  5 on 

the  one  hand  and  3 and  4 on  the  other,  the  lowest  value  of  V_ 

£* 

occurred  during  the  presence  of  the  5-wave  vacillating  pattern. 

The  series  A experiments  were  performed  with  two  different 
cycle  periods  of  T.  (1968  and  984.1  rot).  An  experiment  with  a 
shorter  cycle  period  within  the  same  limits  was  attempted  but 
was  unsuccessful  due  to  the  relatively  large  thermal  mass  of  the 
inner  bath  fluid  and  the  high  heating  and  cooling  rates  required. 
Accordingly,  it  was  decided  to  conduct  an  experiment  with  a re- 
duced amplitude  of  the  imposed  temperature  contrast  variation 
(viz.,  7.8  ^ AT  i 10.0®C)  and  a cycle  period  of  492.1  rot.  The 
results  of  this  experiment  (SE-31;  experiment  B in  Fig.  2 and 
Table  1)  are  shown  in  Fig.  24.  The  experiment  was  begun  with  a 
sustained  period  of  maximum  imposed  temperature  contrast, 
followed  by  36  "annual"  cycles,  after  which  AT  was  sustained  at 
the  mimimiim  imposed  temperature  contrast  for  approximately  2200 
rot  and  then  carried  through  8 additional  cycles  before  termin- 
ating the  experiment.  The  dashed  line  above  the  contour  diagraun 
depicting  the  variation  of  T.  reveals  that,  at  this  high 
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cycling  frequency,  T.  could  not  be  controlled  accurately  enough 
to  make  all  the  cycles  identical.  The  differences  aunong  the 
various  cycles  were  not,  however,  large  enough  to  invalidate  the 
data  from  this  experiment.  The  results  show  that  the  4-wave  ampli- 
tude vacillation,  which  formed  during  the  initial  sustained 
"winter"  condition,  weakened  and  was  replaced,  between  the  third 
and  fourth  cycle,  by  a 5-wave  amplitude  vacillation  which  con- 
tinued for  the  remainder  of  the  experiment.  This  is  not  surpris- 
ing when  we  exaimine  where  the  extremes  of  AT  corresponding  to 
this  cycle  lie  relative  to  Fig.  17  and  recognize,  too,  that, 
when  the  imposed  cycle  frequency  is  large,  the  limits  of  varia- 
tion of  the  internal  temperature  gradient  (5*T/  3R)  are  reduced 
considerably  from  the  magnitude  they  would  attain  in  null  exper- 
iments at  the  two  extremes  of  AT.  It  is  also  worthy  of  note  that 
the  mean  AT  of  this  experiment  is  closer  to  the  4 5 than  to  the 

5 -»■  4 transition  point.  Finally,  we  note  that  in  this  experiment, 
as  in  the  preceding  ones,  is  sustained  at  a lower  magnitude 
during  5-wave  amplitude  vacillation  than  during  4-wave  eunplitude 
vacillation. 

It  will  be  recalled  from  Fig.  3 that  the  series  A and  B 
experiments  were  conducted  in  a region  in  which  the  eddy  thermal 
variance  is  not  necessarily  greater  when  the  imposed  temperature 
difference  is  greater.  The  results  presented  in  Figs.  18  - 24 
reveal  that  changes  in  wave-form,  rather  than  variations  of  the 
imposed  temperature  contrast,  were  responsible  for  the  observed  . 

time  variations  of  the  eddy  thermal  variance.  Fig.  25  (after 


Winston,  Krueger  and  Haines,  1965)  gives  evidence,  however,  that 
in  the  earth's  atmosphere  time  variations  of  the  eddy  thermal 
variance  are  controlled  by  the  seasonal  cycle  of  differential 
heating.  This  may  be  inferred  from  the  figure  because  of  the 
high  time  correlation  between  eddy  available  potential  energy 
and  eddy  thermal  variance  in  the  earth's  atmosphere.  In  order 
to  simulate  this  effect,  we  conducted  two  additional  experiments 
which  are  designated  in  Figs.  2 and  3 and  in  Table  1 as  exper- 
iments C and  D (SE-33  and  SE-22,  respectively) . 

In  experiment  SE-33  the  rotation  rate  was  held  constant  at 
the  same  value  as  in  the  series  A and  B experiments,  but  the 
imposed  temperature  difference  was  cycled  in  the  interval  . 

1.3  5.  AT  i 5.7®C  over  a period  of  984.1  rot.  The  1.3°C  extreme 
is  believed  to  be  close  to  the  marginal  stability  curve  separ- 
ating the  lower  symmetric  regime  from  the  wave  regime  (Fig.  3) . 
Accordingly,  it  was  anticipated  that  would  be  lowest  in  the 
"summer"  phase  of  the  cycle.  We  did  not,  however,  conduct  null 
experiments  or  hysteresis  experiments  to  provide  a guide  to  the 
wave-forms  corresponding  to  different  equilibrium  states  within 
the  range  1.3  ^ AT  i 5.7®C.  Experiment  SE-33  (Fig.  26)  was 
begun  with  AT  sustained  at  5.7®C.  The  fluid  behavior  during 
this  time  interval  was  characterized  by  alternating  mixed-mode 
2uapli‘’:>ade  vacillation,  as  we  might  have  anticipated  from  Fig.  17. 
In  the  type  A experiments  this  was  considered  the  "summer"  con- 
dition. In  the  present  experiment  it  is  the  "winter"  condition, 
since  the  largest  imposed  temperature  contrast  in  this  experiment 
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was  5.7°C.  T.  was  then  cycled  21  1/2  times,  during  which  the 
prevailing  fluid  behavior  was  an  "annually"- modulated  5-wave 
amplitude  vacillation,  interrupted  in  some  cycles  by  6-wave 
or  4-wave  amplitude  vacillation.  As  was  the  case  in  the  type  A 
experiments,  the  appearance  and  length  of  duration  of  wave  forms 
other  than  the  prevailing  one  was  sporadic,  suggesting  once 
again  that  the  internal  non-linear  dyneimics  rather  than  the 
"season"  controlled  the  selection  of  wave  form.  On  the  other 
hand,  the  magnitude  of  V^,  as  well  as  that  of  V^,  was  determined 
mainly  by  the  season  as  it  is  in  the  earth's  atmosphere  (although 
wave- form  influences  are  also  present).  It  is  the  latter  feature, 
in  particular,  which  distinguishes  SE-33  from  the  type  A experi- 
ments. 

Experiment  D (SE-22)  was  conducted  in  the  range  5.5  < AT  < 
10.1°C  with  a cycle  period  of  1320  rot  following  an  initial 
period  of  approximately  1000  rot  during  which  the  imposed  tem- 
perature contrast  was  maintained  at  T = 10.1®C.  The  results 
pertaining  to  this  experiment  are  presented  in  Fig.  27.  Data 
corresponding  to  the  early  part  of  the  experiment  are  missing 
because  of  failures  in  the  data  acquisition  system.  The  mea- 
surements reveal  that  was  larger  in  the  "winter"  phases  of 
the  cycles  than  in  the  "summer"  phases,  as  we  might  have  anti- 
cipated from  Fig.  3.  The  contour  diagram  shows  that  5-wave 
2unplitude  vacillation  occurs  more  frequently  in  "summer"  and 
4-wave  amplitude  vacillation  more  frequently  in  "winter", 
although  there  is  sufficient  intra-"annual"  variability  to 
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classify  the  fluid  behavior  as  indeterminate.  There  is  also  a 
tendency  for  spectral  broadening  during  "winter" . The  only  time 
a given  wave  scale  dominates  through  a full  cycle  is  during  the 
third  "summer".  During  this  period  V„  is  greater  than  normal. 
More  generally,  high  values  of  V_  occur  during  the  presence  of 
4-wave  amplitude  vacillation  and  during  spectral  broadening, 
whereas  low  values  of  Vg  are  found  during  5-wave  amplitude  vacil- 
lation, suggesting  that  the  time  variations  of  Vg  in  this  exper- 
iment depend  not  only  upon  the  magnitude  of  the  imposed  temper- 
ature contrast,  but  also  upon  the  time  sequence  of  wave  scales 
which  develop.  In  this  sense  experiment  D is  a hybrid  between 
the  type  A and  C experiments.  The  fact  that  the  behavior  of 
the  fluid  is  dominated  by  alternate  periods  of  4-  and  5-wave 
amplitude  vacillation  would  not  have  been  anticipated  from  the 
results  of  the  null  experiments  shown  in  Fig.  2.  It  should  be 
noted,  however,  that  the  amplitude  vacillations  in  this  experi- 
ment were  rather  weak,  exhibiting  a variation  of  approximately 
30%  about  the  mean  amplitude  of  the  waves,  in  contrast  with, 
nearly  100%  in  the  other  experiments. 

The  range  of  external  parameters  over  which  the  experiments 
were  performed  was  expanded  further  by  conducting  experiment  E 
(SE-30)  with  a larger  excursion  of  the  imposed  temperature  con- 
trast (viz.,  11.0  5 AT  i 20.2®C).  The  rotation  rate  selected 
was  n = 2.196  sec”^.  According  to  Fig.  3,  Vg  would  not  be 
expected  to  vary  significantly  as  a function  of  AT.  The  results 
corresponding  to  this  experiment,  presented  in  Fig.  28,  confirm 
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this  expectation.  The  experiment  began  with  AT  sustained  at 
11.0°C  for  approximately  1200  rot.  During  this  period  the  fluid 
displayed  a mixed-mode  amplitude  vacillation.  Subsequently, 
when  AT  was  cycled  with  a period  of  2516  rot,  alternate  periods 
of  5-wave,  4-wave,  3-wave  and  mixed-mode  amplitude  vacillation 
appeared.  After  the  fourth  cycle,  some  of  the  thermistors  at 
radius  R = 5/6  gap-width  failed,  so  that  it  was  not  possible  to 
obtain  accurate  measures  of  and  beyond  this  time.  Except 
for  the  general  broadening  of  the  spectrum  in  "summer"  and  nar- 
rowing in  "winter"  there  are  no  notable  seasonal  effects  evident 
in  the  contour  diagram.  Several  winters  are  characterized  by 
4rwave  zunplitude  vacillation  and  several  by  either  3-  or  5-wave 
amplitude  vacillation.  Once  again  we  interpret  this  as  a 
property  which  characterizes  the  flow  as  indeterminate . 

11.  Comments  and  conclusions 

(1)  In  an  effort  to  perform  meaningful  experiments  from 
the  standpoint  of  their  contribution  to  the  study  of  climatic 
variability,  we  have  been  careful  to  make  the  ratio  of  the  im- 
posed cycle  period  to  the  thermal  lag  time  comparable  to  that  in 
the  earth's  atmosphere.  We  have  not,  however,  insured  that  the 
number  of  rotations  per  "annual"  cycle  be  the  seune  as  that  in 
the  atmosphere. 

(2)  Our  null  experiments  have  established  that  the  eddy 
thermal  variance  (Vg)  varies  both  with  wave  scale  and  with  posi- 
tion in  dimensionless-parameter  space.  At  fixed  values  of  the 
imposed  temperature  contrast,  increases  with  distance  from 
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the  marginal  stability  curve.  When  two  different  wave  scales 
are  possible  at  a fixed  point  in  dimensionless  parameter  space, 
the  lower  wave  number  possesses  a larger  thermal  variance  than 
the  higher  one.  As  a consequence  of  these  properties  the  rela- 
tive magnitudes  of  in  the  "summer"  and  "winter"  phases  of 
cyclic  experiments  depend  strongly  upon  the  structure  of  the 
regime  diagram  in  the  region  of  dimensionless  parameter  space 
in  which  the  experiments  are  conducted.  In  the  earth's  atmosphere, 
the  large-scale  eddies  are  observed  to  be  more  vigorous  in  winter 
than  in  sximmer.  Our  results  suggest  that  this  may  be  due  to  the 
particular  location  of  the  atmosphere  in  dimensionless-parameter 
space  due  to  its  special  combination  of  differential  heating  and 
rotation. 

(3)  Our  experiments  have  revealed  that  the  phenomena  to  be 
expected  in  any  region  of  diraensionless-parcuneter  space  depend 
upon  the  direction  of  approach,  the  particular  hysteresis 
properties  in  that  region  of  the  regime  diagram  and  the  rate 
at  which  the  externally  imposed  conditions  are  being  changed. 

Mull  experiments  and  hysteresis  experiments  have  been  found  to 
display  determinate  behavior  at  different  points  in  dimensionless- 
pareimeter  space.  Cyclic  experiments  have  been  found  to  display 
indeterminate  behavior.  There  is  evidence  to  suggest  that  the 
behavior  of  the  fluid  in  cyclic  experiments  tends  toward  that 
in  null  experiments  in  the  limit  of  very  rapid  cycling,  and 
toward  that  in  hysteresis  experiments  in  the  limit  of  very  slow 
cycling.  At  intermediate  cycling  rates  the  fluid  behavior  is 
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no  longer  predictable  solely  as  a function  of  the  externally 
imposed  conditions,  although  some  experiments  displayed  a statis- 
tical bias  for  certain  types  of  flows  in  certain  seasons. 

(4)  The  experiments  in  this  study  were  performed  in  the 
region  of  dimensionless-parameter  space  in  which  the  spectrum  of 
possible  wave  states  is  narrow.  The  range  of  cycle  periods  was 
also  limited  by  the  apparatus  available  at  the  time.  A valuable 
follow-up  to  this  experimental  series  would  be  one  in  which  a 
broad  range  of  cycle  periods  is  employed  in  the  region  of 
geostrophic  turbulence  and  in  the  boundary  region  between  geo- 
strophic  turbulence  and  regular  waves. 

(5)  Referring  to  the  identity 

-Is  = (isothermal  slope)  x 4^ 

d K dZ 


We  may  note  one  significant  difference  between  atmospheric 
behavior  and  the  behavior  of  the  fluid  in  rotating  annulus 
experiments.  In  the  earth's  atmosphere,  time  variations  of  the 
horizontal  temperature  gradient,  and  of  the  related  horizontal 
temperature  variance  are  associated  with  changes  in  the  slope 
of  the  isotherms,  with  little  variation  observed  in  the  space 
mean  vertical  stratification  of  the  atmosphere.  In  the  annulus 
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experiments,  on  the  other  hand,  such  variations  are  accompanied 
by  significant  changes  in  the  vertical  stratification,  with 
much  smaller  changes  observed  in  the  slope  of  the  Isotherms. 
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Catastrophic  Changes  in  Circulation ; 
Flow  Patterns 


Catastrophe  theory  offers  new  insight  into 
fundamental  processes  underlying  laboratory 
simulation  of  abrupt  changes  in  large-scale 
atmospheric  circulation  patterns 


A heightened  awareness  of  the  impact 
of  changes  in  circulation  patterns  of 
the  upper  atmosphere  on  surface  cli- 
mate has  renewed  interest  in  labora- 
tory simulation  of  certain  aspects  of 
these  phenomena.  The  relatively 
simple  gei^lqrsical  flows  generated  in 
a laboratory  annulus,  a ring-shaped 
open  trou{^  have  been  extremely 
useful  in  identifying  and  studying 
some  of  the  basic  mechanisms  of 
these  changes  in  a controlled  setting. 
Such  experiments  are  generally  con- 
ducted by  subjecting  a fluid  within 
the  annulus  to  a series  of  carefully 
controlled  shifts  in  temperature 
contrast  and  rotation  rate.  The  flows 
that  result  cover  a broad  spectrum 
ranging  from  radially  symmetric  to 
turbulenL  These  flows  generally  ex- 
' hibit  fast  motion  in  a narrow  band  of 
the  fluid  when  viewed  (wm  above  and 
are  usually  classified  according  to  the 
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wavelike  pattern  formed  by  this  fast 
stream.  Results  are  commonly  rep- 
resented by  “regime  diagrams”  whose 
coordinates  are  proportional  to  the 
externally  imposed  temperature 
contrast  and  rotation  rate.  The  re- 
gime diagrams  show  the  boundary 
between  ^e  symmetric  flow  and  wave 
flow  “regimes”  and  are  characterized 
by  regions  of  partictilar  wave  forms, 
abrupt  changes  in  wave  form  between 
different  regions,  and  extensive  hys- 
teresis or  delay  in  the  points  of  tran- 
sition between  these  regions.  Thus 
such  diagrams  demonstrate  the  full 
range  of  stationary-state  flow  pat- 
terns of  the  system. 

Laboratory  experiments  designed  to 
simulate  seasonal  variations  require 
a continuous  and  periodic  variation  of 
the  imposed  temperature  contrast 
within  the  regime  diagram.  To  un- 
derstand the  flow  field  which  results 
from  such  conditions  it  is  essential  to 
be  able  to  represent  the  boundaries  of 
the  various  stationary-state  regions, 
the  transitions  between  them,  and 
particularly  the  associated  hysteresis 
behavior. 

Our  aim  here  is  to  present  a new 
mathematical^setting  for  the  discus- 
sion of  these  experiments,  based  on 
Rene  Thom’s  recently  developed 
“catastrophe  theory”  (1975),  an  in- 
novative mathematical  language  for 
the  description  of  abrupt  changes  in 
behavior.  According  ta  this  theory, 
many  discontinuous  phenomena  can 
be  classified  (up  to  arbitrarily  close 
approximation)  as  one  of  a finite  list 
of  standard  behavior  types  called 
“elementary  catastrophes.”  The 
sUrting  assumption  of  the  theory  is 
that  certain  "control”  parameters 
(directly  controllable  or  at  least  ob- 
servable by  the  scientist)  influence 
other  internal  “state”  parameters 


through  some  equilibrium  piinetple 
such  as  the  second  law  of  thennody- 
namics. 

In  our  case,  the  control  parameters 
will  be  the  temperature  contrast  and 
rotation  rate,  and  the  states  will  be 
certain  energies  whose  values  deter- 
mine the  number  of  waves  in  the 
pattern.  Let  us  now  concentrate  on 
the  control  settings  at  which  “be- 
havior” (i.e.  equilibrium  values  of  the 
internal  states)  jumps  discontin- 
uously.  Any  setting  of  controls  de- 
termines a point  in  “control  space,” 
and  those  control  points  at  whic^ 
behavior  jumps  diswntinuously  are 
called  collectively  the  “catastrophe 
set.”  The  theory  then  asserts  that  the 
behavior  of  the  system  near  any  par- 
ticular point  on  the  catastrophe  set 
must  be  like  one  of  the  elementary 
catastrophes  in  the  finite  list.  For  the 
aimulus  experiments,  as  a omse- 
quence  of  there  being  only  two  onttrol 
parameters,  there  are  only  two  ele- 
mentary catastrophes,  the  “fold”  and 
the  “cusp.”  For  three  controls,  there 
are  three  additional  elements^  ca- 
tastrophes, for  a total  of  five;  for  four 
controls,  there  are  seven  elementary 
catastrophes  altogether.  These  con- 
sequences of  the  theory  are  indepen- 
dent of  the  number  of  intenul  state 
parameters. 

The  elemenlaiy’  catastrophes  may  be 
studied  and  anal>'zed  geometric^y. 
With  the  aid  of  schematic  represen- 
tations, the  discontinuous  jumping 
behavior  can  be\‘isualized,  first  in  the 
standard  models  and  then  in  more 
exact  models  of  the  particular  situa- 
tion being  studied.  In  adopting 
Thom’s  theory  for  the  representation 
of  the  abrupt  or  “catastrophic” 
changes  in  flow-field  behavior  and 
associated  hystercMS,  we  hope  to 
provide  new  insight  into  the  funda- 


mental  processes  underlying  such 
evpnis. 

Although  catastrophe  theory  is  used 
here  in  the  context  of  laboratory 
simulations  of  geophysical  flows,  it  is 
by  no  means  restricted  to  such  ap< 
plications.  Indeed,  Thom’s  theory  can 
be  thought  of  as  a prediction  that 
certain  shapes  (catastrophes)  are  the 
most  likely  to  occur  in  nature.  The 
manifestation  may  be  in  space-time: 
Thom’s  original  idea  was  to  classify 
the  forms  most  likely  to  be  observed 
in  early  embryonic  development.  His 
“controls”  were  three  space  coordi- 
nates and  time;  his  “states”  were  all 
the  biochemical  variables  of  the  em- 
bryo. As  we  will  Qlustrata  here,  ca- 
tastrophes may  also  be  discovered  in 
data.  Our  presentation  of  this  new 
application  of  the  theory  is  intended 
to  draw  attention  to  its  potmtial  as  an 
aid  both  in  the  investigation  of  com- 
plex data  and  in  the  design  of  future 
experiments. 

Experimental 

background 

Laboratory  experiments  simulating 
■ certain  aspects  of  large-scale  atmo- 
spheric dynamics  have  been  con- 
ducted by  numerous  investigators 
over  the  lit  twenty-five  years  (Fultz 
et  aL  1958,  1964;  Hide  1958,  1970; 
Kaiser  1970;  Pfeffer  et  aL  1974).  In 
most  of  these  studies  the  annulus 
containing  the  working  fluid,  a liquid, 
was  rotated  about  its  vertical  axis  at 
a constant  rotation  rate,  fi,  while  the 
inner  wall  of  the  annulus  was  main- 


tained at  a cold,  and  the  outer  wall  at 
a warm,  temperature  to  produce  an 
externally  imposed  temperature 
contrast,  aT.  For  different  values  of 
the  control  parameters  Q and  aT  the 
contain^  fluid  will  exhibit  different 
flow  patterns  and  different  relative 
balances  of  fundamental  forces  driv- 
ing the  flow. 

The  streak  photographs  of  Figure  1 
illustrate  some  of  these  flow  configu- 
rations. ’The  symmetric  pattern  of  (A) 
is  observed  at  relatively  low  fi  and 
high  aT.  The  wave  patterns  shown  in 
(B)  and  (C)  are  manifestations  of 
fundamental  instabilities  of  the  flow 
field  and  occur  at  higher  fi.  The  five* 
wave  pattern  will  occur  at  a higher  fi 
with  aT  fixed  or  at  a lower  aT  with 
fi  fixed.  Under  appropriate  condi^ns 
these  wavelike  flows  serve  as  ideali- 
zations of  large-scale  atmospheric 
motion  and  share  certain  of  its  prop- 
erties. The  striking  resemblance  be- 
tween the  idealized  laboratory  wave 
pattern  and  a large^wale  atmospheric 
flow  field  as  represented  by  a weather 
chart  at  a he^t  of 500  mb  b shown  in 
Figure  2. 

A typical  regime  diagram  schemati- 
cally representing  the  results  of  a 
flow-field  experiment  b shown  in 
Figure  3.  The  transition  curves  on  the 
boundaries  between  the  regions  were 
obtained  by  Fultz  et  aL  (1964)  in  an 
experiment  in  which  AT  was  in- 
creased very  slowly  and  fi  was  held 
constant  For  example,  as  aT  crosses 
the  4-3  line  in  Figure  3,  a four-wave 
flow  pattern  breaks  up  and  a three- 


Figurt  1.  Streak  photographs  illuatrato  thre* 
typical  flow  paiums  In  a differentially  heated 
rotating  annulus  of  fluid;  (A)  symmetric,  (1  •• 
1.88  radians  per  see;  (B)  three-wave,  11  ■ 2.19 
radians  per  sec,  (C)  five-wave,  (1  - 287  radians 
per  sec.  All  experiments  were  conducted  with 


a aT  of  20*  C.  In  these  tw»«econd  time  ex- 
posures the  streaks  are  formed  by  aluminum 
flakes  suspended  In  the  interior  of  silicone  fluid 
and  represent  motion  in  the  upper  layers  of  the 
fluid.  The  camera  was  rotated  counts  'ock- 
ivise,  with  the  annulus;  the  flow  and  the  drift 


wave  pattern  forms.  This  process  of 
breaking  up  and  reforming  b a tran- 
sient effect  and  b not  represented  in 
the  regime  dbgram. 

Hysteresis  and 
intransitive  behavior 

The  phenomenon  of  hysteresis,  a lag 
in  effect  when  the  forces  acting  on  a 
body  are  changed,  was  first  observed 
by  Fultz  (1958,  1964;  Kaiser  1970). 
VVhen  experiments  are  conducted  in 
which  aT  b changed  very  slowly, 
hysteresb  b observed  in  the  transi- 
tions betwem  adjacent  wave  patterns. 
Figure  4 Qlustrates  thb  effect  throu{di 
an  overlay  of  two  regime  diagrams 
using  Fultz’s  1964  data.  As  the  over- 
by  dearly  reveals,  there  are  regions  in 
the  repma  diagram  of  a pven  system  - 
in  which  behavior  b bimodah  tiiat  is, 
depending  on  the  hbtoiy  of  the  con- 
trol parameters,  time  are  at  least  two 
stabb  fluid-flow  configurations  for 
the  same  external  conditions.  Such 
behavior  has  been  described  by  Lo- 
renz (1967b)  as  “intransitive.”  Prog- 
ress has  been  made  (Lorenz  1962; 
Cole  1971)  in  attempting  to  exhibit 
thb  phenomenon  in  mathematical 
modeb  based  on  substantially  sim- 
plified governing  equations  of  mo- 
tion. 

We  suggest  an  alternative  method  of 
representing  these  experimental  re- 
sults, using  Thom’s  catastrophe 
theory.  The  model  we  propose  for  the 
qualitative  behavior  of  annulus  ex- 
periments b,  in  a certain  sense,  the 
simplest  possible  model  that  accounts 


of  the  wav*  ptUam  wcf*  site  couoUidorWw*. 
The  radial  lint*  in  th*  photographs  ar*  txltrior 
to  th*  fluid.  (Bxp*rim*t<tal  drtaks  ate  thaaama 
asinPftfferctal.1974.) 
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Figur*  2.  Comparisoa  ahodM  timiUritiM  b*> 
twfndagWirififlOmbpKiiuwpatUniiatht 
upptratinMpbiwof  tl»Nofd»»fB  H>iiiiiphtn 
and  a four*wa«a  pattan  in  tha  laboratoiy. 


for  the  observed  hysteresis  and  in> 
transitivity.  ■ - : ‘ 

The  number  of  waves  of  a flow  pat- 
tern is  a convenient  means  of  class!- 


symmaUie  (low 


wave  Row 


log  Ta  oc  log  n’ 

Figurj  3.  A diagram  indicrl»»af»*.i'>e- 

ary-stata  flow  cliaractaristica  In  an  annulua 
txperimanL  Symmatric  flow,  as  in  Fig.  lA, 
occurs  whan  control  paramatars  ara  to  tha  Itft 
of  the  black  line.  Tha  black  lina  rapraaanta  tha 
boundary  between  tha  symmatric  flow  and 
wave  flow  retimaa,  while  tha  colored  lines  in- 
dicate observ^  boundaries  between  regions  of 
different  dominant  wave  forms  of  tha  flow 
field.  The  number  marked  In  each  r^ion  rep- 
resents the  number  of  waves  apparent  in  the  v 
flow  patte  n,  and  tha  two  numbers  on  tha 
boundarie-  or  transition  curves  between  the 
regions  Iniicate  the  sequence  of  the  transition. 
(After  Fu'.tt  ct  al.  1964.) 


(Laboratory  flow  canditiona  ware  «lma»T  to 
those  in  Fi^  1.  azcapt  Q • 1.95  nsdiboa  par  see) 
In  the  atmoaphera  tha  flow  is  appeoximataly 
pandMtotbaisobaiafthaQowistothacight,  - 


from  high  toka»piaaaura),s»iAspiidinaawiIy 
proportional  to  tha  apactng.  Chaagaa  hi  tl» 
- wava  pattom  have  a aignififant  afSMk  on . 
laigawcalawaatharaadciimataL...-.- 


fying  flowfislds  and  constructing  ta- 
gime  diagrams.  An  intamally  mta- 
sured  paramatar,  howavar,  is  battar 
suited  to  tha  mathematical  modding 
to  be  considered  here.  The  compo- 
nents of  an  atmospheric  energy  cycle 
proposed  by  Lorenz  (1955)  have  bean 
useful  in  studying  long-term  fluid 
flow  both  in  the  abno^thm  (Krueger 
et  al.  1965;  Winn-NieLson  1967;  Oort 
1964)  and  in  the  annulus  (Pfeffer  et 
aL  1974).  Lorenz  breaks  up  the  usable 
energy  into  kinetic  and  available  po- 
tentitd  energy,  each  of  which  is  the 
sum  of  a “zonal”  and  an  “eddy” 
component  The  zonal  available  po- 
tential energy  is  the  amount  of  energy 
which  would  exist  if  -the  masa  field 
were  replaced  by  its  average  along 
latitudinal  circles,  and  the  eddy 
available  is  the  excess  of  the  t'>tal 
available  potential  energy  over  its 
zonal  component  Thus  the  eddy 
available  potential  energy  is  roughly 
that  which  can  be  converted  directly 
into  nonsymmetric  motion. 

The  energetics  of  the  aimulus  exper- 
iments is  not  completely  understood. 
Nevertheless,  some  generalities  are 
accepted.  The  total  energy  of  the  fluid 
(or  any  component  of  total  energy  as 
defin^  by  Lorenz)  is  a continuous 
function  of  time  as  long  as  n and  aT 
are  varied  continuously.  Moreover, 
the  stationary  values  of  some  (per- 
haps all)  of  these  energy  components 
are  correlated  with  wave  number,  and 
there  b a jump  in  the  stationary  value 
when  wave-number  transition  lines 


are  crossed  in  a regime  diagram.' 
These  last  conditions  are  borne  out 
for  the  available  potential  energy 
components  by  ezperiments  tecmiUy  ' 
completed  by  Butyna  et  aL  (in  piepL). 
In  experiments  with  maz^  internal 
theimbtor  probes  (Pfeffer  et  aL  1OT4) 
the  components  of  the  availal^  po- 
tential energy  can  be  meesured.The 
eddy  available  potential  energy  b a 
particularly  sensitive  indicator  of 
changes  in  wave  form  and  hystere- 


Thus  we  can  construct  an  energy  re- 
gime diagram  which  looks  qualita- 
tively like  Figure  4;  the  transition 
lines  now  represent  discontinuous 
jumps  in  energy-stationary  values.  ' 
Our  model  will  offer  a mssns  for  the  . 
repressntaboQ  of  sudx  diagrams. 

• ...  • I • w * V "k 

The  cusp  modfel  - ' - j. 

The  total  energy  (or  any  Lorenz 
component  of  energy)  of  the  fluid  in 
a particular  annulus  experiment  can  . 
depend  only  on  0,  aT,  and  time. 
Letting  E denote  energy,  or  mean 
energy  where  high-frequency  time- 
dependent  behavior  at  foed  aT and 
Q b observed,  it  seems  reasonable 
that  £ should  ^ a smooth  function  of 
these  variables  for  a given  set  of  initial 
conditions.  Thus  it  seems  reasonable 
to  assume  that  theca  u a smooth 
function  /(Q,  aT,  £)  whose  value  et 
a point  in  (Q,  aT,  £)-speca  b the 
tendency  of  £ to  increase  or  deoeese  . 
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“control''  variables  and  S as  a space 
of  "state”  variables.  In  the  applies* 
tions  above,  C u the  (Q,  A70-plane 
and  5 may  consist  of  one  (or  more) 
components  of  the  Lorenz  energy 
cycle.  The  models  we  consider  are  of 
the  form  ‘ 

A«f(c,s)  . 

where  LCXS-^Sisa  continuous 
state-valued  function  of  control  and 
state  variables  and 


/dsi 

ds„\ 

VdT'' 

•'dT/ 

As  used  above,  the  word  “generic” 
applied  to  the  model  means  that 

. f=-VsF 

where  F is  a real-valued  function  on 
C XS  (i.e.  of  control  and  state  vari- 
ables) with  continuous  partial  deriv- 
atives of  all  orders,  and  that  F is 
structurally  stable  in  the  C~-topolo- 
gy.  The  partial  gradient  VsF  is 

/if  li.) 

\dsi ds„/ 

VVe  will  explain  structural  stability 
after  establishing  some  notation. 

For  any  real-valued  F on  C X S,  let 
Mf  be  the  equilibrium  points  of  the 
associated  vector  field,  i.e.  Aff  con- 
sists of  those  points  (c,  s)  in  C X 5 at 
which  VsF  is  zero.  The  projection  r: 
CXS~*C  restricts  to  a mapping  v/-: 


Now,Fis  structurally  stable 
if,  given  a point  p on  Mp,  ^ere  b 
some  < > 0 and  At  such  that  whenever 
G and  an  partials  of  G of  order  are 
within  c of  those  of  F,  then  there  are 
smooth  local  coordi^te  changes  ^ 
around  p in  C X S and  ^ around  «(p) 
in  C such  that  and  ^(Aff)  ~ 

Ma  (near  p).  To  put  it  another  way, 
F is  stable  if  small  perturbations  of  F 
do  not  change  the  qualitative  behav- 
ior of  the  equilibria  of  its  derived 
model 

Theorem  1 (Whitney,  Thom,  Math- 
er). If  the  dimension  of  C is  5 or  less, 
the  structuraUy  stable  functions  on  C 
X 5 form  an  open,  dense  set  in  the 
C'-topology — i.e.  any  smaU  smooth 
perturbation  of  a stable  function  is 
stable  (this  is  the  definition  of  sta- 
bility!), and  any  function  with  k con- 
tinuous derivatives  can  be  C*-ap- 
proximated  by  a stable  function. 

Thus  if  a mo<)el  f of  any  natural  phe- 
nomenon is  given  by  a potential 
function  F,  the  modeler  is  forced  to 
assume  F is  structurally  stable.  To 
begin  wUh,  almost  all  potential 
functions  are  structurally  stable;  and 
having  w unstable  potential  function 
is  equivalent  to  being  unmodelable, 
since  any  error  in  constructing  the 
model,  no  matter  how  minute,  would 
produce  qualitatively  incorrect  be- 
havior. 

Whether  a model  should  be  given  by 
a potential  function  is  another  ques- 
tion. In  many  instances,  such  a po- 


■*V 


tential  function  is  a probabilify  dis- 
tribution which  wants  to  be  mayi-  ,* , 
mized  or  an  energy  function  wlich  • 
wants  to  be  mlxunii^  (Zeeman  1975, 
1976).  In  the  situation  we  are  con-  ■ 
sidering,  the  potential  function  might 
be  thought  of  as  the  “energy  traeopmt 
rate”  function:  the  fluid  in  an  axmulue 
experiment  tries  to  optimize  the 
transport  of  heat  energy  between  the 
outer  and  inner  walls  of  the  azmulus. 

It  is  worth  emphasizing,  however, 
that  the  potenti^  function  need  not 
have  dir^  physical  »ignifif»nw  in 
catastrophe  theory,  its  role  may  be 
simply  to  govern  the  aisoriated  dif- 
ferential equations.  - • 


Theorem  1 gives  the  justification  for 
assuming  that  our  model  b generic. 
The  next  theorem  tells  how  this  as- 
sumptm  leads  to  the  condusions 
reediad  eariier  in  the  pepeA  For  our 
putpoeae  here^  we  can  get  by  with  the 
folloering  concept  of  “aizvalarhy.’* 
Assume  f is  a generic  model  with  as- 
sociated func^n  FH7  X S real 
numbers.  A point  p on  Mp  is  regular 
if.there  is  soms  nmgbborbood  If  of  p 
on  Aff  such  that  -rp  is  one-to-one  on 
If.  CIbe  real  deflnition  is  based  on  the 
fact  that  Aff  is  a differoitiable  man- 
ifold: in  this  context,  a regnlar  point 
of  xf  is  a point  p at  which  tl»  in- 
trin^  derivativeGay  b runsingular.) 
Otherwise,  p is  a singular  point,  or  a 
singularity  of  xf.  Singularitisa  ^ xy 
are  important  for  the  model  .f'- 


because  they  repreeent  pointa  oti  the 
equilibrium  surfoce  where  bifurcation 
is  possible  as  control  parameters  are 
slowly  changed.  Conversely,  if  p is  a 
regular  point  on  Aff,  and  if  If  is  the 
associate  neighborhood,  then  mov- 
ing control  parameters  about  in  x(L0 
will  move  equilibrium  states  about  in 
1/  in  a tottdly  predictable  way,  be- 
cause x|G:l/  -•  x((/)  is  invertible^ 
The  dynamics  of  the  model  keep  rest 
states  on  the  equilibrium  surface 
whenever  possible  Thom  refers  to  the 
imageof  thesingularitin’inder.xas  ' - 
the  “catastrophe  seL”  Thus  in  the 
canonical  cusp  discussed  above,  the 
catastrophe  set  is  the  V-sbaped  curve 
276^  B ^3  jj,  the  control  plweshown 
inFigxireG.  ' 


Theorem  2 (Thom).  If  F is  structur- 
ally stable,  ^en  Aff  is  a manifold  of 
the  same  dimension  as  C and  any 
singularity  of  rp  is,  up  to  smooth 
changes  of  coor^nates,  an  "eleznen- 


I 


£ary  catastrophe.'*  The  number  of 
t,''pes  of  elementary  catastrophes 
depends  only  on  the  dimension  of  C. 
If  the  dimension  of  C is  less  than  5, 
this  number  is  finite,  as  given  here: 


dim  C 

1 

2 

3 

4 

5 

6 

no.  of  elementary 
catastrophes 

1 

2 

S 

7 

11 

00 

Thom  goes  on  to  classify  the  ele* 
mentary  catastrophes  for  dim  C S 4. 
In  particular,  when  dim  0 = 2,  the 
only  two  possible  types  are  folds  and 
cusps.  The  canoni^  cusp  is  given  by 
the  function 


F(a.b.x)°^-a^  + bx 

on  C X S ® l(a,  6)<space]  X |x- 
space]  (o,  b,  z)*space.  The  behavior 
of  the  derived  model 

i = — VsF  = — (x®  — ox  + 6) 


was  described  in  detail  above.  The 
simpler  fold  is  given  by 

F(o,  x)  = Y - ox 

The  reader  may  enjoy  exploring  its 
derived  model.  See  Thom  (1975)  or 
BrScker  (1975)  for  descriptions  of  the 
other  elementary  catastrophes. 

When  S has  many  variables,  say  si, 
. . . , 5n,  the  canonical  cusp  is  given 
by 


G(a,b,si, 


and  the  derived  model  is 

— (si®  — osi  + b) 

h = ~S2 


in  “ -»/i 


i.e.  motion  is  linear  in  all  but  the  first 
variable.  Thus  according  to  Thom’s 
theorem,  no  matter  how  many  “state” 
variables  we  choose  to  observe  in  the 
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annulus  e.xperiments,  it  is  always  at  

least  theoretically  possible  to  change 


coordinates  smoothly  in  and  near  a 
given  bimodal  region  in  such  a way 
that  afterward  one  variable  exhibits 
the  pathological  behavior  of  the 
model  while  the  other  variables  be- 
have linearly.  Some  energy  compo- 
nent or  other  internal  variable  must 
behave  according  to  the  model  based 
on  (Eq.  5);  which  variable  does  so  may 
change  from  cusp  to  cusp.  A fortunate 
property  of  the  annulus  experiments 
is  t^t  one  energy,  the  eddy  available 
potential  energy,  can  represent  the 
entire  family  of  cusps. 

The  regime  diagram 

The  mathematical  theory  we  have 
presented,  together  with  the  experi- 
mental evidence,  leads  to  the  ines- 
capable conclusion  that  the  behavior 
of  an  annulus  experiment  is  closely 
approximated  by  a model  £ » /(Q. 
aT,  E)  where  f “ — V^F,  F is 
structurally  stable,  and  vp  has  a cusp 
over  each  of  the  lines  of  hysteresis  in 
the  control  plane  (regime  diagram),  as 
long  as  ft  and  aT  are  varied  slowly 
and  represent  the  only  control  pa- 
rameters. In  the  context  of  this  theory 
a schematic  representation  of  the 
behavior  observe  in  Figure  4 is 
'shown  in  Figure  7.  Note  that  as  a 
consequence  of  the  theory,  each  hys- 
teresis region  in  the  control  plane 
must  converge  to  a single  point.  Thb 
behavior  is  suggested  but  not  repre- 
sented explicitly  in  the  superimposed 
data  of  Figure  4.  Subsequent  very 
careful  experimental  work  by  Kaiser 
(1970)  has  revealed  that  the  transi- 
tion curves  on  both  increasing  and 
decreasing  aT  do  indeed  emanate 
from  the  same  point  for  a given  hys- 
teresis region.  Additional  discoveries 
by  Kaiser  deal  with  further  details  of 
the  transition  near  the  symmetric 
rt  gion  and  can  easily  be  incorporated 
in  the  surface  represented  in  Figure 
7. 

As  we  proceed  deeper  into  the  wave 
regime  (i.e.  farther  from  the  sym- 
metric regime  boundaries)  the 
structure  becomes  more  complex.  We 
begin  to  encounter  multiple  folding  of 
the  surface,  as  demonstrated  by  the 
sheets  of  wave  numbers  4, 5,  and  6 in 
Figure  7.  We  have  already  seen  how  a 
simple  fold  can  give  rise  to  bimodal 
behavior,  and  multiple  folding  in- 
creases the  number  of  energy-sta- 
tionary values  (or  flow  patterns) 
possible  at  a given  point  in  control 
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Figur*  & Th*  diafiia  ilww*  tfaMT  for  _ 

(Eq.3)forfizMl*.EachpUnra  > constant  is' 
- ODS  of  tho  ptaosa  rsprsaantad  in  Fix.  S>.  Tba 
tuDo^Iow  Haas  art  not  indiostad  hsra,  but  thsy 
arc  tba  saiaa  aa  tbosa  of  Fix.  S,  Motion  of  X is  in 
tba  vertical  diiccfriitowari  the  suTfaf%ascspt. 
over  tbacusp  ration  where  naitioa  is  avwy  bofB 
the  inner  aWt  and  toward  the  outer  sheets. 
The  surface  itself  is  composed  of  equilibrium 
points,withtbainiiarpoidoBofthacusprep> 
iesentinx>instahlaequilibeiaapdtfaarestoftha 
surface  reprssantinx  stable  equilibria.  The 
V-tbaped  curve  shown  in  color  oonatitntaa  tba 
catastrophe  set. 


parameter  space,  as  is  also  obswrved 
in  annulus  experiments  (Hide  and 
Mason  1975).'^en  increased  layer- 
ing is  combined  with  dose  packu^  of 
sheets,  ever-present  perturbations 
may  “bump”  an  energy  state  from  one 
stable  attracting  layer  past  a neigh- 
boring unstable  layer,  at  which  time 
the  dynamics  would  shift  the  energy 
state  to  a different  attracting  layer. 
The  possibility  of  such  “bumping” 
may  give  rise  to  what  we  might  cs^ 
“wave^iumbwvaciilation,”  where  the 
flow  would  be  characterized  by  dif- 
ferent wave  forms  at  different  times 
but  at  the  same  value  in  the  control 
plane. 

As  we  proceed  farther  yet  into  the 
wave  regime  it  will  no  longer  be  pos- 
sible to  represent  the  flow  field  in 
terms  of  single  wave  numbers,  as  in 
the  data  for  Figures  3 and  4.  The 
flow-field  configuration  becomes 
more  complex  and  requires  an  in- 
creasing number  of  modes  for  its  de- 
scription. As  turbulmce  is  ap- 
proached and  encountered  it  becomes 
increasingly  necessary  to  characterize 
the  flow  in  terms  of  a spectral  energy 
distribution.  As  such  complexities 
arise,  it  might  be  necessary  to  con- 
sider Individual  mode  energies  sepa- 
rately or  all  together  as  the  variables 
in  a multidimensional-state  space. 
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Fipuro  7.  A seboiiiatic  •quilibrniai  turfaeo 
reprMont*  tho  •xporimoat^  nnilts  shown  in 
Fit-  4.  Ths  colortd  solid  and  dashtd  linss  out* 
lininn  tht  wavodlow  pattarn  consUtuU  ths 
catastroplia  sat  of  tba  oMdaL  Tht  numbars  on 
tha  surfsca  iodicata  tha  watt  pattarn  obsanrad 
whan  E lias  on  that  portion  (d  tha  surfaoa.  (£ 


Quantitative  modeling 

For  quantitative  purposes  we  need  to 
know  the  function  f with  reasonable 
precision..  Of  course,  one  way  of  ob- 
taining / b to  derive  it  from  basic 
physical  principles.  Another  possi- 
bility is  to  derive  F from  some  basic 
energy-minimizing  principle  or  from 
thermodynamics.  The  latter  possi- 
bility seems  to  warrant  further  study. 
Completion  of  such  a program  would 
yield  a quantitative  model  that  b 
mathematically  and  physically  sound 
in  the  traditional  sense. 

What  we  have  done  thus  far  is  to 
build  a qualitative  model  that  b 
mathematically  sound,  has  clearly 
and  simply  stated  assumptions,  and 
incorporates  well-established  • em- 
pirical evidence.  In  the  absence  of  a 
deduced  quantitative  model,  we  could 
cumb'ne  our  deduced  qualitative 
model  with  experimental  data  and 


logTa  , . _ 

can  rspnMot  th*  a vanabi*  potoatial  Mi- 

•rgy  or  tba  thamut  varianca  of  tha  flow  flakL) 
Behavior  naar  aach  Cttip  b timilar  to  that  of  tha 
aintplar  modal  dbcua^  in  Fig.  6.  Ebpari- 
mantal  raiults  ara  not  accurata  anousb  to  da- 
tarmina  tha  alopa  of  tha  curvaa  ahown  in  color 
as  thay  approa^  tha  cusp  point. 

induce  a quantitative  modeL  (We 
refer  to  such  a procedure  as  “cali- 
bration” of  the  fast  dynamic.)  A 
careful  calibration  might  be  useful  in 
computer  simulations  of  annulus  ex- 
periments, espedally  in  the  “dynamic 
case”  descril^  below. 

Thus  far  our  discussion  has  been 
limited  to  static  modeb:  those  in 
which  control  parameters  are  varied 
very  slowly  so  that  the  system  can 
retain,  or  quickly  return  to,  its  state 
of  equilibrium.  Another  interesting 
situation  is  one  in  which  control  pa- 
rameters are  varied  on  a time  scale 
that  interacts  with  characteristic  in- 
ternal time  scales  of  the  physical 
system.  A problem  of  thb  sort  b 
posed  by  laboratory  simubtion  of 
certain  aspects  of  large-scale  atmo- 
spheric seasonal  and  climatic  vari- 
ability. In  thb  case  aT  b varied 
sinusoidally  at  a fixed  frequency 
while  Q is  kept  constant  Before 


modeling  such  behavior  i t b necessai  y 
to  consider  and  clearly  define  ejl 
characteristic  time  scales  of  the 
problem,  both  in  the  lalxNratory  and 
in  the  proposed  mathematical  settiz^. 
Adding  periodic  dynamics  to  aT  in 
the  static  model  introduces  a second 
time  scale  that  must  be  compared 
with  the  response  time  scale  of  the 
equation  £ = /(Q,  aT,  E).  (In  the 
canonical  cusp^  when  dynamic  conteol 
parameters  are  introduced,  l/e  b a 
measure  of  the  time  reeponse  of  the 
system;  othw  daracterbtic  times  are 
determined  by  the  control  r***t*Tit-  ' 
ters.)  . . 

In  the  laboratory  experiment  there 
are  at  least  three  characteristic  time 
' sca^  that  must  be  coniidaTed:(l)  the 
period  of  the- dominant  physical  ' 
mechanism,  such  as  aT"pl«fa»«lw  ea-  - 
cQlatioii  (Pfrffsr  and  Chiang  1967)  or  ■ 
tilted  trough  vacillation  (FultretaL 
1958);  (2)  “nepoaae”  time  (Lsi.  an 
e-foldiiq(  time  required  to  reettain  - 
equQihrium  after  a change  in  contred 
parameters);  and  (3)  the  external  . 
forc^  period.  For  example^  in  a , 

. particulv  annulus  (Buqma  at  aL,  in  ’ 
prep.)  a chareeteristie  time  in  cate-  . 
goiy  (1),  due  to  amplitude  vacillation, 
b 200  sec.  In  the  same  annulus  the 
thermal  response  time  b approxi-  ' 
mately  1,000  sec.  In  deciding  what 
period  to  assign  to  aT,  we  are  most 
interested  in  numbers  which  yield 
thermal  lags  of  between  0J)85  and 
. 0.17  cycles  of  aT  (one  to  two 
. “months”).  In  thb  example,  forcing 
periods  of  one  to  two  houn  result  in 
lag  times  of  660-860  se^  or  0.18  to 
0.12  of  the  period,  respectively.  Thus 
the  action  (2)  will  be  at  le^  four  . 
times  as  fast  as  (3). 

In  modeling  the  static  case  we  have  . 
filtered  out  the  high-frequency  time 
dependence  (1)  by  coisidcring  mean 
value  of  energy  over  a period  of  va- 
cQbtioiu  These  mean  energies  are  the 
energies  correlated  with  wave  num- 
ber, as  discussed  earlier.  In  the  dy- 
namic case,  it  is  again  these  mean 
energies  which  reflect  the  flow-pat- 
tern behavior.  As  long  as  the  time 
scale  (3)  is  considerably  slower  than 
(2),  it  seems  reasonable  to  build  a 
dynamic  model  by  simply  adding 
periodic  forcing  to  aT  in  the  static 
model.  We  hope  to  carry  out  sudi  a 
program  in  the  future.  (Compare 
Zeeman  19720 

In  thb  connection  it  b interestingto 
note  results  obtained  by  Bur^yna  et  aL 
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' *in  recent  experiments  with  periodic 
variation  of  aT  (in  prep.).  They 
found  that  when  the  ratio  of  thermal 
lag  time  to  cycle  time  is  0.02,  the  be- 
havior can  be  represented  accurately 
by  a static  cusp  model.  However, 
when  this  ratio  is  increased  to  0.12, 
the  flow  begins  to  show  a large  lati- 
tude in  the  point  of  transition  be- 
tween different  flow  patterns  in  the 
hysteresis  region,  and  thus  in  a sense 
the  flow  pattern  appears  “unpre- 
dictable” at  a given  point  with  respect 
to  the  behavior  observed  with  a much 
lower  value  qf  the  time  ratio.  It  ap- 
pears that  above  a certain  value  of  the 
ratio  the  flow  is  able  to  undergo  a 
transition  in  wave  form  anywhere  in 
the  bimodal  region.  In  the  context  of 
the  theory,  the  time-flow  lines  should 
reflect  such  behavior  and  represent 
both  the  history  and  the  probable 
point  of  wave-form  transition. 

Higher  dimensional 
problems 

In  a situation  where  intransitive  be- 
havior is  encountered  but  experi- 
mental evidence  does  not  conform 
with  a simple  cusp  or  fold  model,  yet 
there  is  nevertheless  reason  to  believe 
a generic  model  eusts,  one  is  forced  to 
conclude  that  some  uimoticed  exter- 
nal pa.'ameter  is  operative.  Con- 
versely, we  introduce  another  ex- 
ternal control  parameter,  say  viscos- 
ity, a three-dimensional  control  space 
must  be  considered:  there  is  now  the 
possibility  of  three  additional  singu- 
larity types.  These  new  catastrophes 
will  be  composed  of  families  of  cusps 
together  with  a new  point  at  which 
pathology  is  essentially  three-di- 
mensional. (In  much  the  same  way,  a 
cusp  is  composed  of  a family  of  folds 
together  with  a point  at  which  pa- 
thology is  essentially  two-dimen- 
sional.) 

This  more  complicated  situation  will 
require  careful  study  to  decide  which, 
if  any,  of  the  essentially  three-di- 
mensional catastrophes  is  inherent  in 
the  energ>’  dynamics.  (Since  the  cusp 
is  the  only  essentially  two-dimen- 
sional catastrophe,  this  question  al- 
most answers  itself  in  the  case  of  a 
two-dimensional  control  plane.)  Ob- 
viously, viscosity  is  difficult  to  vary 
continuously  and  slowly  while  main- 
taining constant  (Q,  aT),  and  thuc 
data  pertinent  to  this  three-dimen- 
sional case  will  be  difficult  to  ob- 


In  a situation  where  there  are  com- 
plex data  to  be  analyzed,  Thom’s  ex- 
tremely flexible  theory,  with  its  pre- 
diction that  only  certain  special  types 
of  catastrophes  are  likely  to  occur,  is 
an  invaluable  tool.  In  employing  this 
theory  to  construct  a model  that  ac- 
counts for  the  phenomena  of  hyster- 
esis and  intransitivity  in  circulation 
flow  patterns,  we  have  demonstrated 
only  one  of  its  possible  uses.  Thom’s 
theory  is  widely  applicable,  and  al- 
lows new  insight  in  constructing  glo- 
bal dynamical  models  from  experi- 
ment^ data  as  well  as  pointing  the 
way  to  the  design  of  more  effective 
experiments. 
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Abstract 


-^Rotating  annulus  experiments  are  described  in  which  the 
imposed  thermal  forcing  is  varied  periodically  by  means  of 
cyclic  variations  of  the  cold  inner-bath  temperature.  In  the 
way  of  background,  the  results  of  conventional  annulus  experiments, 
in  which  the  imposed  temperature  contrast  is  not  varied,  and  of 
hysteresis  experiments,  in  which  the  imposed  temperature  contrast 
is  varied  slowly  enough  to  be  considered  quasi-static,  are 
reviewed It  is  found  that,  in  the  region  of  dimensionless- 
parcuneter  space  in  which  our  experiments  were  conducted,  the 
flow  configuration  in  conventional  and  hysteresis  experiments  is 
determinate  in  the  sense  that  it  depends  only  upon  the  point  in 
dimensionless-parameter  space  at  which  the  experiment  is  conducted 
and  the  direction  from  which  this  point  is  approached.  The  flow 
configuration  in  experiments  in  which  the  imposed  temperature 
contrast  is  varied  in  a cyclic  fashion  is  found  to  be  indetermin- 
ate in  the  sense  that  it  changes  from  one  state  to  another 
sporadically  and  cannot  be  predicted  solely  as  a function  of 
iseasont . 


Cyclic  Variations  of  the  Imposed  Temperature 
Contrast  in  a Thermally  Driven  Rotating  Annulus  of  Fluid^ 
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